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INTRODUCTION 

Each subroutine in the Scientific Subroutine Package is documented by means of 
a flowchart. In the Guide to the flowcharts each flowchart is listed under the 
general heading "Statistics", "Mathematics" and under a subheading relating to 
its basic function. 

SYMBOLS USED 

Figure 1 illustrates the various blocks used in the flowcharts and their particular 
meaning. Lines connecting these blocks are made up of asterisks or periods. 
Arrows showing the direction of flow are represented by a "V" for downward, a 
left parenthesis "(" for right to left, and a right parenthesis ")" for left to right. 
In some cases an X is used to show flow direction. Special symbols are used 
for mathematical operators within the flowchart boxes. These are: 

* for multiplication 

& or + for addition 

— for subtraction 

BLOCK NUMBERING 

Numbers appearing inside the connector circles refer to the chart block location 
on the page; these are numbered vertically zero to nine and horizontally A to C. 
Numbers shown outside the blocks correspond to FORTRAN statement numbers 
in the prc^rams. 

EQUATIONS 

Many of the flowcharts contain references to equation numbers. These equations 
can be found in the mathematical description for the particular subroutines in 
the User's Manual (H20-0205). 



( ) 



Enter or exit block 



Processor block 



Modification block 




Decision block 
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Call to subroutine block 



Connector to block within page 



Connector to block on another page 



Figure 1. Flowchart blocks 
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JKCA) 


K 


* 


♦ 


liO* 




IhAN 


Vh-X 




« 


« 


* 



* fij $v «««« *««<-«* 



*«**«««*««««**« 



« KLt-'LACt VHAX**««I* bO« 

• av A( iJI ♦ * « 

* * «** 



rtl flKST • ♦♦* 

tfAlilAbLE • * * 

titr^tNT IN *«*«!* C0« 
Uu-iLKVATiL'M ♦ • « 

MATklX * ««* 



SUOKQUTINt XA^iZ 



* jti) •*♦)♦ 



» * 

tALCUL«ilt * 
K^LATlVt « 

IrOUATlCM ^1* 



« « 

♦ LALCULAlt ♦ 

• MtANS * 

« * 



* i.ALi:ULATt * 

* st.inOaiu: * 

« UcViATlUtt^ « 

* (tt^UATiON 41* 



*•»*♦****«"••*** 



«««4«««««»«4»«* 



SUttROUTINE SUBHX 



ENTER ; 



«*« 




♦ 




* 




*«* 


• • 


* 








•NO 


• • 


* eo«i* 




LAST 


COLUMN 


•!• A3* 


* * 


* 


UF 


MAJHIK A 


• 


• • 


•*« 




• 




• 




♦•« 



• ««* •^•i**** •• « 



* SET L'C LL'O 



«**•««•*•**« «4* 



* • 
•COUNT N0HZ6H0* 

# COOES IN « 

* VECTOR S « 

• * 
•***«*«««*«*4«* 



FOtt EACH 

CCLUHN OF 

DATA HATRfX 

A, DO 

FOLLOhING 



*«*•« **** ****** 

EXIT 
*************** 



*************** 
*** * * 

* « * 

* A3****l* SET I«l 

* • • 



*** * 

* * * 

* ^^s****^* 

* * * 
*** 



*************** 

* • 

♦ ♦ 

ADO 1 TO L * 

* 
*************** 



c 



*************** 



• AOO 1 TO LL • 



*************** 



**** *********** 

* * 

* HOVE DATA * 

* FRON AILl TO* 

* 0(LLI * 

* • 



*i********************* 



•vts* * 

IS I EQUAl TO •)* BO* 
f LAST ROW • * * 
* * *** 



*************** 



AOO t TU 1 ****l* A4* 



*************** 
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SUOROUTIME COft^e 



ENTCit 



•«•«••*«*«•••«• 



* • « Oasi«VATIONS« 



• CtJ«**«l« SUeTKACT I • 

• * * FAUN KK « 



MMICIH6 
SrORACE 



• ottfivATIaN 






MC OAT* 



• CM.CU(.«TE • 

• NUNS. SEE • 

• CaiMTIOH Oi* 



•««*«««x» 



« 

• afyfiRriotis • 

• food IWANS • 

• « 



MIC 



S Of • 



AuSfSSS? ev 

fSMTION III 



****«*•«**•«*««•**•*«•< 



SUM OEVIHION* 

Him • 
TMmwMr • 



«••«• 



• SUM CROSS * 

• rSwUCTS OF * 

• OCVIATIONS • 



**♦♦•♦♦*•*•••*• 



* CALCULATi 

• HK-N-K 



♦«* •OATA * 

tS*«M|* GET AN • 

• * ossSitrAriuN * 



««*«**•**««•««• 



* SUM THE • 

• OtSEKVATION * 



• CALCULATE * 



«**««**«******4 



2US V 

• CALC SUMS OF* 

• CROSS * 

• PRODUCTS Of • 

♦♦*•♦•••••»♦••• 



• CALC CORK. • 

• COtfF. IN • 

• STMHETRIC • 

• MATRIX R. EO* 

• Ul » 



* « 

* CALCULATE * 

* STAMOARO * 



•COrr DIAGONAL* 

"ciiSi": 

iSU**«***« 



IJT » 

•DATA 

• mSMVATIl 

* 



• MSEIIVATION • 



EXIT 

♦•♦»>»—♦ » <♦*** 



IWSHVATIONS 



.SUN OJJIAT 
* TEMPUUK 




AIlONS *•••)• CO* 

— I TfNf. • * • 

MEANS * **• 



tl 



SUBRUUT[Nt HISK 



*t ************* 

ENTEK 
*************** 



*************** 

* * 

* SUM CROSS ♦ 
•I* PRODUCT Of * 

♦I-TH AND J-TH* 

* VARIABLES • 
***********^**^ 



************ 



** 



• COflPUTC ♦ ♦ *, 

*)*$KEHNESS AND ••♦♦♦•♦♦•♦•**00*l 

♦ KURTOSIS. ♦ ♦ * 

* EQ. (11-121 * ♦* 



******•«. 



*************** 



*************** 



*«•**«*•«**•**♦ 



NON-HISSING 
OATA GREATER 
» THAN 2 * 



********t 



* *oo 1 1 

* ♦ 



« * * 

« * * 



,* **************** 



* * 
«CCHPUTE MEANS* 

* FOR * 

* VARIABLES. * 

* EQ-Il) • 
*************** 



*************** 



* ADO 1 TO K 



*«*•***•******* 



************ 



* lER-l 

• * 
« * 
•♦♦♦*♦*••••♦**♦ 



**********.* 



*************** 
EXIT 



,************** 



.****•*••****•***«. 



*************** 

* * 

* SET J=l 

* 

* ♦ 
*************** 



*************** 

* COMPUTE SUM * 

* Of CROSS • 

* PRODUCTS OF ♦ 

* DEVIATIONS. * 

* EQ.(2) • 
*************** 



ARE eOTH 
VARIANCES GT. 

THAN 
* 10**(-201 * 



********^ 



*************** 
♦* • ♦ 

•A4***}* SET 1=1 

* « * 
** * ♦ 



*************** 
** * * 

* « * 
•A5***l* SET K=l 

* • * 

*************** 



** * *YES ••• 

• * • IS OATA ♦ ♦ ♦ 
•A6«l XtKjII •! Bl* 

• • • HISSING ♦ • * 
•• * * ♦*♦ 



* * 
*» CI* 

* * 



• IS 

NON-HISSING 

♦ OATA GREATER 
♦ THAN 2 • 



*«*«****. 



*************** 

* 

* IER"l 

• * 
*************** 



*** 
* * * 
*«•«) DO* 



* * 

♦ ♦ 

♦ ieR= 

* * 
*************** 



* ♦ * 

••••I C6* 

* * * 

*** 



,*****************t 



*************** 

* * 
*COMPUTE REG. • 

♦ COEf. AND * 

• INTERCEPTS. • 

* eg. 14,5,7, «)• 
*************** 



*************** 

* * 

* COMPUTE * 

* CORRELATION • 
*COEFFlCIEttTS.* 

* EQ.OI * 
*************** 



c 



^***************** , 



♦YES •*• 
IS DATA ♦ ♦ • 

• *• 



1SSIN6 



IS VARIANCE 

GREATER THAN 

10»*l-20l 



**•**••*. 



*************** 

* • 
♦COMPUTE STD. ♦ 
♦ERROR OF REG.* 
♦COCF, EO.IA • 

• AND 9} * 
*************** 



*^ ***************, 



*************** 

* * 

♦SUH I-TH ANC * 

• J-TH • 

• VARIABLES * 

• * 
*************** 



••••*«*•****•*♦ 



*** 



* ** 

♦* 



♦ ieR=z ****} 00* "ce*} 
« «****» 

* * *** * ** 
*************** * 



♦YES ♦♦• 
* * * 

*) DO* 

♦ «•* 



.*«***♦*♦*♦♦«*♦«**, 



«♦♦«*«*•♦♦♦♦*♦* 

♦ ♦ ♦«♦ 

♦ SOH SQUARES ♦ ♦ ♦ 

♦ OF I-TM ANO ♦♦♦♦) BO^ 

♦ J-TH ♦ ♦ ♦ 

♦ VARIABLES ♦ **• 
♦♦♦♦♦*♦♦♦♦♦**♦♦ 



♦♦•♦♦♦*♦•**♦♦♦• 

♦ ♦ ♦♦* 
•COMPUTE STO. • ♦ ♦ 

♦ DEVIATIONS. **♦*) CO* 

• 60.U0) ♦ * ♦ 

* * *** 
*************** 



*************** 

* * 

* * 

* AOQ 1 TO 1 

* • 

* • 
*************** 



*** 

* * * 
****) A5* 

• * « 
*♦* 
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:»UtlROUTINt OKOtR 



CNTtB 

0* ****** ******* 



*************** 

♦ COPY CORK ♦ 

• COtFf 4X VS * 

* Y) IN ftY<JI * 

• FtJK RET- * 

* VWUBLES • 
0**0 *********** 



•*4« *********** 

• COPY • 

* INTER-CORR * 

• COEFF OF • 
•RETAINED INO.* 

* VMIABLES * 
*************** 



« S(ORE * 

« OeP€NOENT * 

* VARIABLE * 

* NuneEK IN * 

* ISAVECK*!) ♦ 



*************** 
EXIT 
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SUURUUTINC miLTK 



ENTFM 



*•« • • 

♦ • • CALCULATE^ ♦ 

• * «L4U4T1UN U2)* 

•*♦ * ♦ 



* CALCULATE 
KEI 

see 

* eOUATIONU) « 



• « 

« CALCULATE • 
•STDf FlUtUR OF* 
«£STlHATb. SeE« 
*EtiUATIUN (10l« 



* CALCULATE • 

* COEFF Of * 
VOETMMINATIOM* 

* . SE» 4 
« MUATION 4AI« 



• * 

• CALCULATE • 
« F-VALUL. • 

• EQUATION lai* 



««4*4«««***« ••* 



♦ CC6FF, __- 

* EQUATION 121* 



• • 
*&tOHE RESULTS* 

• OF * 

• CALCULATiON • 

• IN ANS(l-iOI« 



••** •««4«*** ••« 

• ♦ 
« CALCULATE * 

• INTEKCEFT. • 
« SEE EaUATICH* 

• (3) * 
•••«**•• ******* 



V 
EXIT 



( 



• CM.CIH.«IC • 

• «3|MTIoN 161* 

• • 



• « 



••MM********* 
. • 

• C<IC1N.«TE • 



• « 

* VUti MCC OF * 
.•eSTIHATE. $Ef« 
I* (OIMTION 111* 

•««««****«««4*« 



•«*«**•***«•••« 




♦***l* 44* 



M 



SUBROUTINI: CU<I< 



• CfMCKATE « • 

• VMI«S(.CS. • 

• MIMTIOH III* 

• * 



• CM.CUI.AIE * 

• Mf*NS. • 

• CeiMTION 121* 

• • 






• CM.CW.*;! • 

• COUELMIOII • 
•COtfFICIENlS.J 



* OEVIATlOMS. • 



.aj.«i- 



•••**••*«*«•«• 



6X1 T 
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SUBROUTINE STPRG 



fcNTEK 



•* 



«**•**•**•**•«* 



COHPUTE P, * 

•BO»*»I*PI>OPORT10N OF* 

* * • SUH OF SO. * 

*• ♦ REDUCED • 



♦♦ ♦ ♦ 

* • *COMPUTE REG. • 
•CD***l»COeFFlCIENTS * 

* * *SEE EQ. 1101 * 

«* « * 

*************** 



*************** 

• ♦ 

♦ IMITIALIZt • 
♦CONSTANT AND ♦ 

* OUTPUT ♦ 
♦VECTOR. lER^O* 
*************** 



!•; 


P GREATER 


* 


* 


• 








•1 


(.9* 




INPUT 


* 


* 


* 


• 


CONSTANT * 






k** 



****•««****•*** 

* COMPUTE THE * 
♦INTERCEPT SEE* 

* EQ, Till * 

* * 



*•*»*****•***** 

* FIND « 

* DEPENDENT • 

* VARIABLE* « 

* FORCED AND ♦ 
♦OELfcTED VARI.* 
**«««*••*••*«** 



*************** 

* STORE * 

* VARIABLE • 
•SELECTED AND * 

* SUM OF SO. • 

* REDUCED ♦ 
*************** 



**•***♦♦*♦♦**♦♦ 

* * 
•COMPUTE STD, • 

* ERRORS OF * 

* REG, COEFF- • 
•SEE EQ. 1121 « 
••*****«****«*• 



***>*•*******•* 
« • 

* FIND THE • 

* MAXIMUM * 

♦ NUMBER OF • 

• STEPS ♦ 

•««****•«*****• 



*«•*******«***« 

* GET * 

* CUMULATIVE ♦ 

• SUM OF SQ. ♦ 

♦ REDUCEO SEE * 

• EQ. (3> ♦ 
*******•*•***•* 



•*•******•***** 

• * 

• COHPUTE T * 

• VALUES SEE ♦ 

• £0. tl3) ♦ 

• * 
*************** 



ARE ALE 

FORCED 

VARIABLES 

ENTERED 



*«•****••«*•*** 

• GET * 

• CUMULATIVE « 

• PROPORTION • 

• REDUCED SEE • 

• EQ. 141 * 



**««********«** 

• PERFORM * 
•REDUCTION TO * 

• ELIMINATE * 
•VARIABLE EQ. • 

• (14-161 • 



4* ««*•** ******* 

* SELECT fHE ♦ 

•NEXT VARIABLE* 

• TC ENTER ♦*♦♦•*••••, 
♦AMONG FORCEC ♦ ♦ 

♦ VARI. * • 
«*««*******••** • 



**•••**«••**•* 
* 

♦ COMPUTE 

* MULTIPLE R 

• SEE EQ. <5) 
* 
*****•••*•*«♦• 



c. 



*************** 

• • 

• COMPUTE ♦ 

• ADJUSTED • 

• MULTIPLE • 

• CORRELATION • 



, ****«***•«****«*•, 



««*«****««*«*** 

♦ SELECT Th£ ♦ 
•NEXT VARIABLE* 

♦ TO ENTER ♦ 

♦ AMONG * 

♦ NON-FORCEO • 
*******•«••**** 



♦••♦•♦••*♦•♦*♦• 

• COMPUTE F ♦ 

• VALUE FOR * 
•ANOVA SEE EO.* 

• 16) * 
••*•••*•••***** 



«*•***«««•***«• 

• • 

• COMPUTE ♦ 
•ADJUSTED STO.^ 

• ERROR OF ♦ 
'• ESTIMATE • 

*•*«***«*««••*• 



*(*♦***«***«**** 



is sum of sq. 

reduced gt. 

Than o 



*«««•«,«*. 



•««««***«****** 

* • 

♦COMPUTE STO. * 

• ERROR OF * 
•ESTIMATE SEE • 

♦ EQ, (71 • 

****«***«•••**• 



*••****•••••••* 

• STOUT • 

• CALL OUTPUT • 

• SUBROUTINE ♦ 

• • 
***♦******♦**♦♦ 



* DOES 

REDUCTIUN 

• EXCEED TOTAL 
♦SUM OF SO. 



♦NO ♦♦• 
• * * 

♦ ) 80^ 

• • 
• •* 



•«*«*•««•**•«*• 

• INVERT ♦ 

* 01 AGONAL • 

• ELEMENT Of ♦ 

♦INPUT MTRIX * 

* SEE EO. (6) * 
*************** 



• •NO •♦• 

• HAS ANALYSIS • • ♦ 
TERHINATEO IN •! A4^ 

♦ STOUT • ♦ • 

* * ••• 



*(****««********♦* 



**«**•*•****•** 



♦*♦ 
• * * 
♦♦•♦I C9* 



*•*******«**«** 



*«•*«*****•**«« 

• * •** 

• DIVIDE THE • * * 

• ROW B¥ THE ••••) CO^ 
♦PIVOT SEE EQ.* • • 

• (91 ♦ ♦♦♦ 



«* «****«***••***• 
• • • 

•C9^l^ EXIT 
« • * 

*♦ *************** 
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SUBRQUTIN€ PROBT 



ENTtR 



* * * 
*B0***1* SET IER»3 

* * * « 



* * • 
•*«*! D8* 



♦COHPUTC l-TH * 
I* WORK I NO * 

* PROBIT SEE * 

• EQ. (6-7) * 



IS I»L*5T 
SII*<ULUS 



*•***•**, 



••••«••*••***** 



••••«««••«***•« 



IS ANY 

OBJECTS 

RESPOND GT. 

*THAN TeStEO* 



****••**, 



*************** 

• * 

• SUN * 
*INTERNEOIArE * 

* RESUtTS * 

* * 



************* 

♦ 

* AOO I TO I 

* * 



* * * 
«*«*) CT* 



,««••**•«*•«*««« 



IS NO. OF 

OaS£ LEVELS 

6T 2 



*************** 
♦ IER«4 



« * • 
••••) 08* 
* « * 



IS I=LAST 

STIMULUS 

LEVEL 



**««*«**, 



* C(WPUTE AO. * 

♦ SUM OF • 



.**«••*«**••«««***, 



• • * • 

• SET IE««l ♦*♦•! DB» 

• * * * 



20 V 

«•«««*« *4«***«« 

♦ COMPUTE * 

• PRGPOKTIONS • 

* Of OBJECTS • 

* RESPONDING • 

• SEE EQ.dJ • 



***«»«*«*•*«*•• 



• ADO 1 TO 



tm**** ********* 



.***••••• 



«**•«*««. 



.****•***•**••*•*♦. 



[S THERE 

HlhUS DOSE 

LEVEL 



•••*«***. 



* * 

* CONPUTE * 

* PROaiT REG. * 

* COEFFICIENT • 

* SEE EO. 1 21 * 



•♦»•••**♦*♦••*♦ 

♦ CONPUTE NEW * 

• ESTIHATE OF * 

* PROBIT REC. * 

* COCff- EQ. ♦ 

• (111 * 
*m*******mm**** 



m****** ******** 

• « *** 

• * ♦ « 
•STORE SXX IN ***•) 68* 

• OD • * ♦ 
» * •** 

M************** 



*t***************^ 

♦ ♦ 

* • 
80* ♦ 



-*«***************. 



* * 

* * 

* SET 1ER=2 



** 



••»*•*•••••■ 



*** 

* * * ** 
***•} 08* ** 

• « • ** 

• **« ** 

* ** 
** 



*************** 

• • 

• COMPUTE THE • 

• VALUE OF • 
•INTERCEPT SEE* 

• EQ. 13) * 
*************** 



*•••*******««•* 

* ♦ 

• COMPUTE NEW * 

• ESTIMATE OF * 
•INTERCEPT SEE* 

♦ EQ. (12 1 • 
*•**•«***««*«•* 



*«***•«*« ****** 

* * 

* STORE THE * 
♦INTERCEPT AND* 

* REG. COEfF. * 

* * 
*************** 



_ •**••******«•**«*, 



*************** 

* COMPUTE * 

* EXPECTED ♦ 

* VALUES OF ♦ 

* PROBIT SEE * 

* EO. (41 * 
***•«***«•••*«* 



*****•**•*«*«** 

« * 

• SET 1*1 

* « 
********+***«*« 



**«*******•*••* 

* * 

* COMPUTE CHI * 

* SQUARE SEE * 

* EO. (151 * 

* * 
*♦*♦•••♦•«••♦** 



*•«••********** 

* * 
* 

♦ SET 00=0 
* 

• ♦ 
****«****«*•*** 



*************** 
** • COMPUTE NEW * 

* * * VALUE OF * 
*C7***)* EXPECTED • 

• * • PROBIT SEE ♦ 
*• • EO- (131 ♦ 

**•«****•**•*•• 



****«*****«**•« 
♦FIND OEGRf" 

• OF FREEOp 

• FOR CHI 

• SQUARE Ste 

• EO. <lfc| 

«****«****•*«*• 



GREES • 

eopM * 

CHi ♦ 



******••««*••*« 

« * 

* CONVERT * 

* STIMULUS * 

* LEVELS TO • 

* COMMON LOG • 

•*•*••••*•**•*• 



*•«****•***•••• 



*************** 



•««i*«««******** 

• FIND 0, * 

* DIFFERENCE • 

* &ETUEEN OLO * 

♦ AND NEW • 

• PROBIT ♦ 
****** ********* 



«« *****«•*******« 

* * * 
•08*)* EXIT 

** 4*4**«*«*««***« 



*VES ♦♦♦ 

* • * 
•I Bl* 

* * 
**• 

*** 

• * 
♦I ♦ BO* 



***•****•«*•*** 



*COMPUTE I-TI 

• WEIGHTING 

• COEFF. SEE 

• EQ, (51 * 
*******««****•• 



*** 

* * * 
**♦*) CO* 

• * * 
*** 



*************** 

• * •** 
*SUM SQUARED Q* « * 

♦ IN SXX SEE *♦**! DO* 

* EQ- (141 * * » 

* * *** 
«***«»•«***•*«• 
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^UttKUUriflt CANGK 



tNTfK 



« pMTinoh • 

* CuiaftCLAIICN » 

• MATRIX • 
« EtiOATIOh Ml* 



• * 1-ti IJ IS 


*NU 


ttU*)* NUM6EK lif 


«*«••*«•*••« 


* * CRITCttlAl 


• 


••• • 


* 



exir 



*«««•••*«•«••*«*««**•«' 



• CALCULATE * 

* EICCNVALuESx* 
« ElCENVECTUftS* 



*NRC]CT • 



«•***«••««•«••« 



••4«**«** *«***« 



* AOU 1 Tli I •*•♦)• AS* 



««««*•«*••••«•« 



sn I - I 



n 



• • 

• C«.CULArE * 
CMONICAL * 

CljllltEI.*TICIl « 

• EOIMTIOM <SI« 



• CM-SOtMRE * 

• EgtMTlcin i«>* 



«*•«««•*««•**•• 



**«*•••*••••• 



»»»» >*#»»*♦**** 



• * 

• C«LCU1*T( • 

• cacrricitKis* 

• EOIMTICHI lai* 



• • ««• 

• CALCULATE • • • 

« LFFT lum *•*•■• »a* 

• ClJEFFICIfNIS* • • 
« EOUATION 191* **0 
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^Ufli<UUTINI AWOAI 



6NTFI* 



«•**««•«******* 



«••«•«••••«••*• 



••*•««»•••••«•« 



* * * 



« • 

• HOWE OATA T0« 
•UPPHt PART Of* 

• ARftAV X « 

• • 

***«****4**«*«* 



* * 

* « 



« CAtC STEPS * 
*USCO TO STOKf* 
* DATA. Sf£ • 



***«*****•**••* 



* AOO i TO « 

* i(4JUNT(Jl • 

* * 



• IMITIAiUC • 

• KOUNT4JI-1 • 

• FOO DATA « 
« SUBSCRIPTS * 

• fOft J-1,K • 



* « 

* • 



««*««**•««•« •** 



«***«*«*••***«* 



*4««««4««**«««« 



AOU 1 TU J 



**«*«••******•* 



• • • 

« « « 



•*•*«•«•*•««*«*•*•** 



*•• « • 

* « •CALC LOCATICK* 

* A6«***l* ILt TO STOKE* 

* • •A DATA. SEE * 
•«« « EQUATION Ul* 

**««•«*««•«• • *« 



«***«*****••*** 



«*««***«*•**««* 



««« 
* * « 

« * « 



•«•***«**«*****«**«* 



• * 

• REOISTftieUTe* 

• I-TH DATA • 
«POIMT TO X(LI* 

• « 



**«««**«****♦*« 



400 L TO I ••••I* A6» 



•*•******•**«** 



•••*•«*****« *4* 








*NQ * 


• 


-0 


•1* 


Ol* 




* * 


« 


• 


«*« 


♦ 






• 


*** 


• TCS 


* 


* 


**•••••**«« 1* 


dO* 



19 



^U^KOUVINk AVCiL 



CNTFR 



***« 



•Yt-S » » » ♦NU 

J=N« ♦♦««♦♦»♦♦«♦♦♦ « CU*1* FIi'^iiH ALL ***♦•♦♦•♦•« 

• » ♦ * * t^ACTUK ^UNS * 



*************** 
» ♦ 

* C«.CULATt • 

* CuMTf^OL F0« ♦ 
♦cNO Of FACTUrt* 

* kUN * 



* * * * 

* * *** 
*************** 



tXlT 



*««*««**«*«*«**«***««*« 



««****«******■•*«*•««*** 



♦ • 

♦ LASTSl ll=L+J* 

♦ **HERfc L 15 • 

♦ LOCATION OF * 

♦ LAST DATA • 
*************** 



* 

* ♦ 

*************** 



* * *** 

* * * * 

* ADO I TJ I **♦*!* At,* 

* * * * 

* * 9** 
*************** 



*************** 

♦ * 

♦ LASTSU) = ♦ 

♦ LASTSlI-ll ♦* 

♦ ISTtfU). fCK* 

♦ I«2.K ♦ 



*** 



*************** 



* CmCCUCAIE * 

* „i****i* ^N*X(LL1-SUM* 

* * * \U0 STGAf iT* 
*♦* * IN X<LLI • 

*************** 



167 

*************** 
*** * CALC OATA ♦ 

* * • RANGe * 

* C3****t* CONTSOL. C =* 

* ♦ * L - LAST&in* 
♦♦* * + ISTEPUl * 

*************** 



*************** 

♦ * 

• IS TCP I II IS * 

♦ I-TH STtP ♦ 
*USfcO TO sruiiF* 

* DATA * 
*************** 



*************** 



*^* ************ 



C GP.tArtK ♦**•**♦♦•♦»♦♦ 

THAN U * 



( 



* 

• SET 1=1 

♦ ♦ 
*************** 



*************** 



* • •SET L=l. * 

* A6****>* LL"l. SUH=C, * 

* ♦ • NN=LEVEL<I) • 
*** * • 



« * * * 

* ADD 1 TO J ♦*♦*»* A3* 

* * * * 

* * *•* 
*************** 



********************** 



*************** 



♦REPLACt L AMD****}* A7* 
• LL bV C+l ♦ ♦ ♦ 



*************** 



************** ***4l***** 



*************** 

♦ * *** 

♦ AOU ISTEPU »♦♦*♦>* AT* 

♦ TU L AUO LL • • * 
« * ««* 
*************** 



*************** 
*** * * 

* * * 

* /il****}* SET J>1 

* • * 

•*♦ ♦ ♦ 



* 


• 




167 




♦ 




*** 






•NO 


♦ * 


L- 


LASTiUI 


♦ 


• C3* 






* 


• « 




* 




«** 


* 


* 








* * 




**♦ 




******* 




* « 




**** 


• CO* 








* * 








«*• 



*************** 



* A««*««1«A00 DATA kILJ* 

* * • TO SUH • 

*** , * * 



*•*•«•****•«** 



• AGO ISTEPUl****)* 6C* 

♦ TO t ♦ * ♦ 

* * *«* 
*************** 



i. 



20 



SUflftuUriht MEAMO 



ENTER 



**« • « 

• * • CAtl, flKST • 

♦ t)a*«*«)* AND iCCUNO ♦ 

* * * JiVIStikS * 
«*« * * 



PERFORM t4€AN 

SQUARE 

OPERATION AS 

FOLLOWS 



SECOND 
OIVISOK IS 
I He bEGk^ES 
oe FREeUOH 



« * 

• INITIALIZE 

«waRK AAEA AND 

• CONTROL FOR 

• LAST DATA • 



♦ nivmt i-TH ♦ 

* ^H ftV fiKsr* 

♦ OIVISOK • 

• « 



* SET NN * I 

« 

« • 



Trie AttllVE 

rttSULT IS SUN 

OF SQUARES 



•** « * 

« • *CALC POSITION* 

• A*****|* IN nUTPUT * 

• • * AREA TO SUM * 
•♦• ♦ NN-TH DATA • 



• « 
« OIVIUE 1-TH * 
«>UN BV SECOND* 

• UlVIS'Jk * 

• • 



« * 

* SQUARE NN-IH* 

* DATA AND * 
^ACCUMULATE IK* 
•THE ABOVE LGC* 
*««*•••««*•**•« 



THE ABUVE 
RESULT IS 
HEAN SUUAKES 



*«*«******«***4 

* * **• 

* « « ♦ 

* ADD 1 TO «N ••••1* A4* 

* * * * 

* * *** 



3UT^^* n« am 
OUTPUT IS 2 
TO THE K-TH 

PGUER 



**4«**«*4**«**««««»***i 



EXIT 

****«**•* ****** 



*«**••*««•«**••**«*••** 



•******«****«** 



* AUO 1 TO I **•*^* H-J* 



•««««««••«««««* 



•*»•••***«•«••« 



«*« 
« * • 
♦I* ttO* 

* « * 
•** 



21 



^UekOUTlNt UHATM 



ENTER 



•• • • 

« * CALC POOtEO • 

• ♦ hjEtkix* • 
* eouation <3)« 



•*« 



« IttlllALIZE 

« UORKINC 

* STOIUCC 

* * 



EXIT 

••♦♦♦•♦*••••••• 



«««*4««***«4««* 



* SET K - 1 



«« ••«**«*««*4** 



••« * SUM K'TH • 

* « GKQUP DATA • 
A3****1*AN0 CALCULATE* 

• • MfAMS • 



•• 


« 




SUftTRACT • 




MEANS FRCH * 




DATA * 




♦ 



( 



•«*«•«♦•«««««** 



• CALC sun Of « 
•CKOSS PCOOUC1* 



• » 

» ACCunuitTe * 

• SUM OF CROSS* 

* rRODUCT or * 

* otviArioMS • 



• ' LAST 

CKour 



« » » » 

* AOO I TO K •***)* A3* 

• • • • 

U^MU******* 



a. 



22 



SUeftOUTINE DISCR 



ENTER 



*************** 
**• *FOR EACH CASE* 

• * 'IN EACH GROUP* 

• B0****1* DO fHt • 

• • • FULLUWING • 
*•* * SIEPi • 

*«•****«**«*«*« 



• « 

• CALCULATE * 

• COMHON N£ANS* 
•fOR VMIABLES* 

• EQUATION I 11* 



* * 

* CALCULATE * 

* JISCRIHINANT* 

* FUNCTION ♦ 

* EQUATION 151* 



* CALCULATE • 

* GENERAL * 
♦MAMALANOeiS 0* 

* SQUARE • • 

* E0.I2I * 
«««•«*•«*•««•* * 



•««••*«****«*** 



••«*«**«**4«*4* 

• • 

• CALCULATE * 

• DECREES OF « 

• mEEOOH * 

• * 
*«*««**«****•** 



« CALCULATE * 
•PROBALITV FOR* 
* LARGEST • 



*•*•«• ******* «• 

• * 
« 

• SET K - I 

« 

• ♦ 
«**««*•*««** *** 



«*«*******•**«* 
EXIT 



•«*•*******«*«* 
•4* * CALC * 

• * • COEFFICIENTS* 

• A5**«*l« OF 0ISCR1HI-* 

• • •NANT FUNCTION* 
•«* * E0.I31 * 

««*«*«*•«•*««•• 



« * 

* CALCULATE * 

* CONSTANT * 

* EQUATION (4>« 



*•• 

•ves* ♦ 

•)• BO* 



«****«•**«**••* 



• ADO 1 TO K M**!* AS* 



•*••••••***••«« 



23 



suaoouTlNi; th/icc 



ENTER 



••»««*«••«•««•• 



«••**«****•*•*« 



EIG£N««IUE 
LESS TIMII 
CONStltNT * 



• • 

* C<1CUL«T€ * 

•rckHEJlTtiE oe; 
•K £tCEN««I.UES* 

«♦**♦♦♦•••••••• 



*•««•«*•« A************' 



c 



l**********««* 



/too 1 TO « • 



CIVIOE I-TH 
ElitHVALUE ■>: 



• • • 

MO I 10 I •«•»♦ «2* 



24 



SU6IIOUTIN6 LOAC 



fcNT» 



SCT 1-i 






•••••• 



NULTirtT l-tM» 
nCENVKTOII * 

■T M • 



V 

exit 



*m*m*»*** * *******»***** 



MO 1 TO I •»*•• *» 



25 



SUBRUUTINE VArtHX 






♦♦* ♦ * 

* » ♦ • 

* e;»**»*»* ADD I TO NC • 

* ♦ * ♦ 

*«* * * 



♦ ♦SECTION (I.-3I* 

**« « Of eUUAflQN * 



* * 

* str €PJ( = 

V0.00L16* LL = 

* K-1. NV = 1 

* * 



« ••♦ 

NC *»EATtR •>♦ Cei* 

THAN 3 * * • 

* «»* 



»*♦ ♦CALCULATfc C-ii* 

* * IHtlA, •il'* » 

C;4**«)« TheTA, i^^ • 

• ** * mi-l Hi 



***•■****»» 



• * 

♦ SET NC = O, 
« CONS « 

♦ 0-7071066 

• ♦ 



« * * 

« ^2****}* SET J = 1 

♦ ♦ * 



*««* 



«**«««* 



♦CALCULAIf tOb* 

♦HHI, SIN ?HI.« 
♦Sfcf tOOAIlbNS* 

♦♦•»*♦♦*♦♦«♦*♦* 



♦ • 

* CALCULATF • 
•COHHUNALlTieS* 

• SEE EQUATION* 

♦ {2} * 



*************** 

*** * * 

* ♦ * 

* Bl****!* SET Kl = J*l 
« * « 



* ♦ 

• PLKfurIM • 
*KUTAIIOK, Stt* 

♦ tOUATlONS * 
*t26> ANU i2?l» 



• CAtCUlJkTE • 

• NORMALIZED * 

• FACTOR * 

• MTRIX. SEE * 

• EQUATION 131* 



m*0*4i* ********* 

*** * • 

« * ♦CALCULATE NUH* 

* e^4*««)* AND OEN. SEE* 

* • * EQUATIUN <7)« 
*** * • 

♦♦♦♦*»*•••♦•♦•• 



#** * * 

« « * « 

* C4****)* ftDU I TU Kl • 

* * * * 
**♦ * ♦ 

*************** 



c 



#«******•*«*•*« 



• <S****I« HDD 1 TO NV 



HUH ttSS THAN .»»»»»•••••»• 
DEN * * 



*** 

♦ NLi • * 



««««••«•**•**•* 



•••••**♦♦••*•*• 

• CU-CULATE • 

• VAHKNCE FOR* 

• FACTOn ♦ 

• MATRIX. SEE • 
» Ca. (41 ♦ 



*¥E5 
NUH CREArER .«••♦«*«» 
THAN OEN • * 



♦**♦♦****♦*•••♦ 
« • 

« * 

• AOO I TO J • 

♦ * 
» • 
♦**••♦♦♦♦***•** 



« * A30 
« « **• 

»»ES» • 

HV EREATER OR •)• C8» 

EQUAL TO 50 • '^,* 



• • 250 
« * ♦*• 

« *N0 • • 

NUH ♦ OEN •)• C0« 
• LESS THAN EPX • *.„* 

, , ... 



* *** 
*NU * * 
J . K-l .1. 83. 



• . 

• PERFORM • 

• CONVERSENCE • 

• TEST. SEE * 

• EQUATION <5|. 



A * 

• rERFORH * 
.OPERATIONS IM....|< 
.SECTION IB-11* > 
. OF EQUATION . 
*************** 



t** 



$$t 



430 

* OENORNALIZe * 

* ca***«l* EACTOK * 

♦ • ♦ HATKIX. SEE ♦ 
• ««. *EQUATiaN 126)* 



****«*«««•««•«««**** 





* 


*** 

•NO * * 


.tIIPIeo 




*i* 82* 




m * * 




• 


*** 


• 






• • 




«** 


• YES 




• « 








* • 






♦** 



***** **** ****** 



* REREORH 
♦OPERATIONS It. 
•SECTION (8-21* 

• OF EQUATION • 
«••••••*•*•«*•* 



26 



f. ci» 



. CALCULATE » •••• 

. FINAL • * '^ 

»C0MHUNAL1TIES»*»«I» AO * 
* SEE EQUATION* . 2* 

. (29) • •*** 

....*..*....•♦♦ 



SUhftUUTlNE VARHX 



m**m *********** 
**** * CALC • 

• * *OIFFERENCE IN* 

• 40 *«*|*COMMJMALITIES* 

• • * S£F EOUATICH* 
**** * (301 * 

*************** 



#« ^^^0*^******* 



27 



SUBROUTIME AUTO 



( 



EHTEII • 

• 



• CALCULATE • 

• AvetASE OF * 
•TINE SE*IES A* 



*•«««> 



* CALCULATE • 

tASifl??!!*!?:.: 

• • 
****««* ***** *** 



EXIT 

*************** 



( 



28 



SUMOUriME MOSS 






29 



SUeROUTIW SMO 



ENTER 



•«*•«««•*«*•*«« 



•*«**«•« «•** **« 



MOOTHED 
* SCRIES R • 



•*«4 **«***•**«* 



*««•••«•**•«••* 



c 



30 



suexouTiNE Exsm 



ENTER 



• ARE •YES 
COEFFICIENTS •«»»•»»«•• 

* A« B« C CIVEN « 



* •• c - 

• EO.ISI-III • 

« « 



FOLLOUIHe FOU 
ALL IMT*. I • 



• rORECAST ONE« 

• fEIIOO MEAD* 

• Ea.iii • 



•WIMTE COCFF.* 
*USINS cumcHT* 

• «. i. c • 



31 



SbdROUTlME KOLHO 






* • 

*ltR<i tRKJk eiAb* 



•«•« 






•««« 




• • 






• • 




• A3 • 






• A4 ♦ 




« « 






m « 




•«** 






«••* 




X 


>• 




X 




• 


• 


• 


« TRAN&FORH X, 




« 


CALL ikHifW 


• 


• CALC. PIXl 




* 


FOR PAOB *> 


* 


« 




* 




• 








« 


* 










I *«*• 




* 




• 


• 








..X* G3 


• 




, 




• 


m 








««• 


■ 








CAUCHY .*. 






X 




Bi *, 




«««««{!<,«»••«•*•«« 


• • •• 




• 






.* FIRSI « 


. VES 


* 






.« 


* 


P-l-P 




*. Z£RQ .• 




* 






• • ■* 




« 






*. .« 






• NU 






• 





• RtOKOtR X»S ♦ 

♦ rtLlNUTOHlC • 

* iNCiiEAStNi^ • 

* • 



ircKATE X 

• cALCULAriONS ♦ 

• THRJUCH ULOCK * 
...X*27 ARC FOR EACH* 

• XUI * 

• * 



» • 

♦ * 

* CALC. P<X) • 

« • 

♦ * 
..X* G3 * 

* • 



u3 *. 

*'* tUtifQlHTS 



««»*C4*««****«* 

• * 
« EX[T • 

• « 



.0U6 



♦ » 

, NO K • • 

,* ...X«S£7 ERROR FLAG * 

X • * 

- • • 

• * 

• IK • 

• • 
• ••« 



••♦♦♦£ ^•••••••♦♦* 

IcALC. SAMPLE OF* 

• FOA POINTS ♦ 

• 1-1,1 • 
« ♦ 





!...!■ 


...cos H2 




....z. 


...005 K2 




... .3. 


...ooa ea 




.***'•* 


...C05 03 




....^> 


...005 FJ 




• « 




♦•♦♦♦ 

•00 5* 
• H2* 

• * 
• 


• 0^ •x.. 
« « . 

♦••* , YES 

.* •• 



• K2* 



«*•«• J 2«« ******** 



* * 

* ii^ *x.. 

* * . 

«««« . YES 

.♦. 

KZ *. 

.* *. 



• * 
« * 

• CALL PIXl *. 

• * 

• * 
•«•••««•********« 

•005* 

• FJ ♦... 

• * . 
«•** 

USER A 

•«v*«F3********** 

• • 

• CALC PCXl Ifl ♦ 

• USER UtFlMEO * 

• PDF * 

• ♦ 



c 



* CUHPARE PtXl- * 
•SAMPLE OF HUM * 

* PRfcVIOUS ONE •K. 

* KEEP • 

* j;4=HAK.0lFF. * 
««•**•««•**«««**• 



* 

* 
• *«« 



* * 
« • 

* * 

* * 
*«*«««•«•*««***** 



,X* 02 * 

« * 



it X 

««««•« 3 •«*«*««««« 



««*« *««••«*«•««*« 



• A4 • 

* * 



c 



32 



&UdAOUTIK£ KOklU 



• tNTfcK K0LH2 * 

• * 



«AcOKOEA X*S AflO* 

• V«S, MONOTONE • 

• INC. * 



• FOR ALL X,Y • 

• CALC OIFf'. • 
«a£TH£Elt SAMPLE • 

• O.F.'StO * 

• « 



• K£TA|N MAX 

* UlFf=« 

« 



•t 3 *«•*««***' 



«*««***♦•***•••*' 



•4f 3 •••*«•* 



««*•••••*•* 



P-l-P 



« * 

• EXIT • 

* • 



33 



SU8A0UTINE SNl RN 



» * 

* ENTER SMIftN * 

* • 



* « 

* * 

* SMIRN - * 

* * 



03 *. 

« 

X LI 1 



^*mm*i:3** ******** 

♦ * 
♦L*XI= SHIRN FOR* 

• )C LT I SEE * 

* MATH 0£SCR. • 

* * 



♦ ♦ 

♦ SHIRN =1 • 

• * 

• * 



♦L*XI= SHIRN FOR* 

♦ X NOT LT * 

• I SEE MATH • 

• OESCR. • 

♦ • 
***m ************* 



«***J3*««****** 
• EXIT • 



( 



34 



SUBROUTINE CHISO 



*•*•••*«*•***** 



** «**« •*****«** 



* BO**** I* S£T ERROR * 

* ♦ • coot - 2 * 
*** * * 



**** **••«•«**** 

* INITIALIZE * 
« ERROR CODE 

* AND 
*CHI-SOUARE TO 

* 2ER0 * 



80 V 

••«**••*****«** 
*** • COftPUTE * 

• * TOTALS OF * 
Bl****>* COLUHHS AMO * 

* * ROMS <Eg, * 
**• ♦ 2-*l • 



* • 

* COMPUTE * 

* DEGREES OF * 

* FREEDOM. « 
« lEOUATION II* 
••*••«****•**** 



**«**«****«***« 

• * 
•COMPUTE GRAND* 

* (£0u8(thN 41* 



•NO 

DECREES Of •*****♦♦♦**•♦ 

FREEDOM - ♦ • 



•YES 
2 BY 2 **•****••***« 

SPECIAL CASE • 



^m************* 



****«****•*•*«* 

• COMPUTE • 

• CHI-SaUARE * 

• FOR OTHER * 

• CASES. * 

• (EQUATION 61* 



EXIT 

«« «••******•••« 



*••**»*«******* 



• ««****** ****•« 



*««4>« *•****•* *•••«*•*•* 



*********************** 



*************** 



••*«•*•««*•**«« 
^fXIT 

•*«*•*••••*••«* 



*••••*****«****•««*•**: 



120 V 

•**••••*•••••*« 

* COMPUTE * 

• CHI-SOUARE * 

• FOR 2 BY 2 « 
« TABLE. * 

* lEOUATIOM^SI* 
*•*•**•*••— ^^^ 



•**•*•••***«•** 

EXIT 
•*«**•••*•***** 



• * 80 

• * **• 
1/5 OF CELLS •NO • * 
CONTAIN LESS *)* Bl* 

THAN FIVE • * • 

* • •«* 

• • 

« * «** 

* VES * ♦ 

*********** I* ao* 



35 



SUttROUIIIfC KKANK 



EMTES 



* eu***«l* SET L ■ 1« 

* « « > N«L 



•♦* •Tit • 

• • mm**** •******•* 

• CO**'**)* CQMPUlfc TA • 

• • • ItUUAIION II* 

• «« « * 



DATA AAHKfcO 



«••«••«««*♦•••* 
*«ANK * 

••«••««•***•••* 



*««*•«••«••«*•« 



•• •« *m*»******* 

• RJMK OAT* «••♦♦••♦•♦••« 

* VECTOR B * 

« • 



««*««4««««*«*******' 



• MOVE RAMEO « 
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♦ tC.iiJ * 



If LPTUfc - 
5, CALC 
StCuM., 



»««««««•**««*« 



UPTll^h = i ♦•***♦*.*•♦*•* 



•KCH£Ki tt,»^J» ***** 



««4«.«»*«v««*«** 



V 



♦NO * • 



♦♦♦♦♦*♦♦••••••• 



*««««««*««*«••* 
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SUeRUUTINE Trt^si 



€NrCK 



IbU 

*♦• ♦ • 

• * • CALC T-VALUe* 

• BU»***I* *N0 OCCfttbS ♦*••♦*♦♦♦•♦••••• 

• • * Uf FKEConH ♦ 

♦♦♦ • tg, (bl-t?) • 



««•«•*•«••••«** 



**«*•«*••**«•*• 



I7d 
*•* « * 

♦ * ♦ CALC T-VALUE* 

• Bi*»**»* »N0 OEGRECS *•*•••*••• 

* » ♦OF FREtOOM ♦ 
••* * tw- 181-491 ♦ 



««*««•«**«**•*« 



•••****♦♦♦•**•♦ 



•**••••*•«*♦**« 
«** • * 

* * • CALC r-vALue* 
« Bd*»**i* AND DEGREFS • 

♦ ♦ ♦Of FftEEOOH * 



«( •»««•**•«•••«**•**«*• 



OPnOh *¥FS ♦ * ♦ 

GREATER THAN ♦**••♦♦*«♦*••♦• 63*«>* EXIT 

A ♦ • ♦ • 



« 


• 180 




♦ ••• 




•VES* • 


OPTION 


= < •I* b2* 



c 



« • 2UU 

* • •«* 

•VPS* ♦ 

OPTION lESS •!♦ Hi* 

THAN 1 ♦ ♦ • 

• « *•* 



•*«*«4««**«*«** 

♦ * 
« CALCULATE * 
•VARIANCES EO.* 
« (2J * 

• * 



OPTION » I ♦♦«•*♦•••*••• 



• « •*« 

• CALC T-VALUE* * • 

• AND DECREES •♦••)* 63* 

• ( OF FREEDCN * * * 

• eg. 131-1^1 • ••• 



«««*«««««*«*•********«* 



* • 170 

* *«* 
•NU ♦ • 

OPTION ' 2 *1* Bl* 

• « • 

• *•« 

« * 

* « «*« 
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SUBRGUTlNe eiSER 



TO HIGHER 
•CATEGORV OF tt* 



•««•••*«••«••*• 



•NOTRI * 

•~~ COMPUTe • 

• OnOINATE OF • 

• HOW • 



• • 
*COHPUrE MEAN • 
•OF COKTJNUaUS* 

• VARIABLE A • 

• SEE EQ. Ul • 



* • 

* COMPUTE R. • 

* BISERIAL • 

* CORRELATION • 

* SEE EQ. (31 * 



* • 

•CCHPUTE SIO. • 
•OEVlATlCh OF • 
•A SEE EO. IZi* 

» • 



•COMPOTE STO- • 

• ERROR OF R • 

• SEE EQ.. 1*1 • 

• • 
••••••••«»••••• 



« FINO Pt * 
•PROPORTION OF* 
« CASES IN * 



«•••«*••••«•*•• 



• CONPUTEO • 

• RESULTS IN • 
•OUTPUT AREAS • 
•*««•«•••***••« 



«««*4«*«**«*«*« 



•FIND g« O-l-P* 



•••••••*4****** 



• * • 

* • « 



••*•«•*•******• 



••••#•••• 



*•«•••*••«•••*• 



• • « 

• •*• 



,***•••••*•••**«•*, 



53 



SUBROUTINE ^Hl 



•CONSTRUCT A 2* 

• X 2 • 

• CONTINGENCY • 



IS ANY CCLL 
EQUAL TO 



*•**«***. 



•••••••«•**•**• 



«•*•«**«•' 



••*•*•••**•*••* 



, ♦•••••••••••♦•♦' 



« COMPUTE rHI • 

• COEFFICIENT • 

• SEE EQ. Ill * 

• • 
******** *•••*** 



••***••***••••* 

• * 

• COMPUTE CHI • 

• • 
•**•••***♦««•• • 



«««*••**•••**** 



KAXINAL PHl 

* COEFFICIENT • 

• • 
****•••*•«•••*« 






•(*«*•***•••*•*•*. 



*******•••••*«« 

EXIT 
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SU6ftOUTIN€ POINT 



A************** 

£NTER 



** 



«***»**«*****•• 



♦COMPUTC HCAN • 

*BO***l* OF A BELONG * 

• • ♦ TO LOWER • 

*• *CATECOAV OF B* 



* * 

• IER=0 * 

• * 

* « 
*************** 



*44************ 

♦ • 
•COHPUTE POINT* 

* 6ISER1AL • 

• CORRELATION * 

• SEE EQ. U) • 



• * 
*CCHPUTE MEAN * 
•OF CONTINUOUS* 

• VARIABLE A * 

• SEE £Q. Ul ♦ 



• * 

• COMPUTE T ♦ 
*RATIO SEE EQ.* 

• U) ♦ 

• ♦ 



«♦•♦•*••«••♦♦♦• 

• * 

•CCHPUTE SID, • 
^DEVIATION OF ♦ 
♦A SEE EQ. 12)* 

• * 



* • 
•FIND DEGREES * 

* Of FREEOOH ♦ 

♦ SEE EQ. T5J • 

• * 



• FIND NO, OF * 

• CASES IN * 
•HIGHERi LONER* 
•CATEGORIES OF* 

• B ♦ 
•*•**••*••••••* 



•••^•••••••••*« 

EXIT 

*••••***«••••«* 



* ^VES 

IS ANY CASE * 
IN LOWER CAT. 
OF B 



**««••«•. 



••••••*••• *4*«* 

« 

* I£R«l 

* • 
« • 

**••••••*•*•••* 



*•• 

• • • 
*••*) B4* 

• « * 



,♦•••*♦♦•••*••••♦«. 



• •NO 

* ARE ALL CASES * 

IN LOWER CAT- •••*••••, 

• OF B • 



••••♦••♦♦•••••♦ 



««*«***••«•**•• 



,*««**««•«••«***•«. 



• * •** 

•COMPUTE MEAN 4 * * 

• CF A ttELONG ••••! 80* 

• TO HIGHER « • * 
•CATEGORY OF B* *** 

••••«•«•••****• 
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SUBROUTINE TETR* 



•CCNSTRUCT * 2* 






• continSem 

• TAM.! ^o 

• INPUT 0«T 



0«T* » 



•**•*««•«••«*«• 



• SET EURO* 

• CODE, IE>a 



cowyTE n * 

flONS • 



• COWUTE r 

• TOTALS ES 



* COMPUTE P2 * 

PROPoifiaNS 

"OH TOTALS 



••*•••«••••*•*• 



IS THERE A 

StTISFACTORr 

ANSWER 



•••«****«**4«*« 



•SET IE-} OR ** 

•»V SUeROUT I NE •••*♦•»*♦•* 

• POLRT, • 

• »-RS-10*»75 • 



.••*••••«*•«••*•• 



S^ANoHl 
ERROR 



•*•«•****** 



5Sl.lS:i»5 : 



•(•••«••••*••«**< 






•NCTRt * 



«••«••«••«••««• 



*NOTIll • 

* • 

* OftOIMATC OF • 

* m.c, *I P2 • 

* • 






AltC THCME 
LflfLC 
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SUUOUTINE S«>T£ 



ENTER 



••«•*«•*«••*«•« 



•PROPORTION Of* 
*OVING SEE Ea.> 
• ill • 



» • 

• INITIALIZE 



••••*•**•••*••« 



• conruTE • 

•PROMRTION or* 
•SURVIVING SEEP 

• ES. I» • 



MP************ 



•PPPPPPPPPPPPPP 

• P 
P COMPUTE P 
PSURVIVtL RATEP 
P SEC EQ. 1*1 P 

♦ • 
PPPPPPPPPPPPPPP 



PPPPPPPPPPPPPPP 



« •••••••••«•**•**•. 



• i'jric'i-tffiT • ppppppp. 

P « M'DRRMN P 



PPPPPPPPPPPPPPP 



IE*2 



PPPP) tTP 



PPPPPPPPPPPPPPP 



IS iPlNnjT 
VALUE OF K 



PPPPPPPPPPPPPPP 

p P 

p P 

P ADO 1 TO I ;• 

p p 

p p 
PPPPPPPPPPPPPPP 



.ppppppppppppppppp. 



.ppppppppppppppppp. 



PYES ^PP^ . 
ARE DATA * P P P 

"'siSWlt pppppppp.PtTPlP 

no Efl.^ . ••„• • 



P ARE D^ 

p Accini 



EXIT 

PPPPPPPPPPPPPPP 



pppppppp< 



PPPPI (T* • 



PPPPPPPPPPPPPPP 



P| MP 
p p p 

Ppp 
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SUBROUTINF HCPY 



$f ************** 
FNIFR 



«♦♦*♦*«♦♦***♦♦* 



***t*******i*** 



H:** *********** 

* * 

* COPY MATRIX * 

* MFHENT !.1Y * 

* (.IfMFNT * 

^mt ************ 



*************** 



*************** 



o 
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SUBROUTINE RCPY 



*«»*«******«**• 



***•*««»*•** 



• » 

• INITIALIZE 

* ?OINTrR AT 
♦FIXST CLEMENT 

* Of kOU * 



««* * ♦ 

♦ * • GET fLtHENT • 

* AZ****»* FROM INPUT • 
« * • MATRIX * 

•»* ♦ ♦ 

♦♦****•*••****♦ 



• MOVE fLfMENT* 

• TO OUTPUT ♦ 

* VttCTOR ♦ 

* * 



LAST FLFMFNT •*«♦»♦*•**♦** 
Of ROH * * 



• STEP POrHTER* ♦ ♦ 

• TO NEXT •***!♦ A2* 

• ELEMENT OF ♦ * * 

• ROU ♦ *•♦ 



♦•*•*♦**♦*****»•*•*♦••• 



EXIT 
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SUBROUTINE CCPV 



CKTfR 



• • 

• IMITIALIZE 

• POINTtR *T 
•FIRST ELEMENT 

« OF COtUHN • 



*•• • « 

• * * GET ELEMENT « 
« A2****l* FROM INPUT • 

• • • HATRIX • 
««• • • 



*«******••****• 



* MOVE ELtHENT* 

* TO OUTPUT • 

• VtCTOR * 

• * 



♦ves 

tAST ELEMENT ♦♦•♦♦•«♦•*♦♦» 
OF COLUMN * 



c 



• « «•* 

• STEP POINTER* • ♦ 

• TO NEXT •••*!• A2* 

• ELEMENT OF • * • 

• COLUMN • *** 
•♦•»♦*•♦♦*•♦•♦* 



*«««««««««***«*******•* 



EXIT 

*««4«***«4««««* 
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SUFKnUTINf OCPV 



CMTfl! 



« INITUllZE 

• POIHTM- *1 

♦rir.si rifMEMi 

* ♦ 



••«v*««4«*«*«** 



» • • OCT 01»G<»*1» 

• A7*»*l« (LtHENT HKIM 
« ♦ • INPUT MTRIX* 

*** # * 



*««**•*•*«*•«** 



» TC OUTI>UT » 
• VICTIM • 

» » 



•YtS 
L«ST ILIMCHT .««•.»♦••»*•« 

nF oiActiNii » ; 



W 
*♦*«*****♦*♦••* 
« * *** 

* STeP PUINTER* * ♦ 

* TO NFKT ••»*)• A?» 
« fLFHEHf OF • • » 
» OIAGtiNAL « **• 
9***9********** 



*********************** 



EXIT 
*************** 
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_^^8^^- 



SUBROUTINE XCPV 



****«««**«*«*** 



«««««« «*«**•« 



««**•*••«••*•» 



• INITIALIZE 
*ARR«V POINTCB 

• AT L.tC IN 

• ftATftlX A • 



•*«««*••«*•««•• 



« A3***«l* G£T ELEMENT • 



*•«««•««•*«««•« 



**•*•******««** 



••*«•«••««««««* 






•NO 

LAST ELEIKN1 •*«»»•••»••». 
IL2iK2l • 



EXIT 



m99*******************' 



««««•«•««*««••« 



STEP PaiNTER* 
TO MEXT 
HEMENT 



••♦•I* A3* 



•••*«*«•••••*•* 
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SUSRCUriNb HSTA 



ENTtK 



« « 

• INITiALIZf 
«ARP.«V PCilNT£B 

• AT FIRST 

• CLfthT * 



«**•«**«•*«*•** 



* A2****t* GET ELEH6NT * 



44** *«#*«****«« 



*YtS 
UUTPUT HATftlX ♦•♦*•••*•••*• 
CtNEKAL * ♦ 



tLFNEhT iK 
LOWFR 
IklANOLE 



•I •«*««««««««***«****** 



DFf-OIAG 

ELENEMT OF 

CI AGONAL 

UUTPUT MATRIX 



* • 

* FORM ELEMENT* 

* OF UUTPUT • 

* HATKIX • 
« * 



* i **** ♦***•♦♦♦♦••♦••••• 



EXIT 

I 



«4««4««««* «««««* »•««««« 



* * •** 

* STEP POINTER* * * 

* TO NEXT ♦♦**|* A7# 

* ELEMENT * « « 

* • •*« 
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SU!iK'}(ITINl' l.fC. 



ENTI K 



■IS (ess THAN *»*J*t«:»««»« 



* « 

* r,VNf:KAL • 

* HATP 1 X • 



«««««««««•*«««•* ::*«*** 



• * «*« 

* CIWPUTE ♦ * * 

♦ •iUuSCRiPT F(ta****I« A..* 
■« SVMMhTSlC * * * 
» KAIHIX ♦ f** 



C 



««•«»««* *ir «*«*«« ««***?* 



*♦*♦•♦♦****»»♦• 

* * 

* cnwPuTf: • 

*SUftSr,RIPT FOrt» 
« OIAGfWAl. * 

* HATRIX • 



*1 ******** *«*«***« 



•*« ***•***««*«**** 

« « * 

* A&««1* *^XIT 

m * * 

*** *«««*«**«*•**** 
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SUSROUTIME CONVT 



ENT 



IS HOOE'l ♦«*••♦•♦••♦•* 



* CONVERr * 

* NUMBERS FROH« 

« S1N(U.E ••♦••*«•••••• 

« PRECISION TO* 

* DOUBLE • 



♦•♦*••«•**•••*♦*•*•* 



♦ CONVERT • 
« NUmeRS FftOH« 

* DOUBLE * 

• P^tECISION TC* 

• SINGLE • 



*| *««««*««***«*«•*** 



exiT 
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SUiSfCIITlNF ftF'<*'' 






• iS Hont = 1 **♦***»**«*** 



LONVtRT FRUH* 
SINCLf TO *♦********•*« 

OIHENSIUN ♦ ♦ 



««««***»«****** 



♦ * 

♦ CONVERT FRt»M* 

♦ UnUBLE TO * 

9 simuE • 

« DIMENSION * 
*«*♦♦«•**♦♦*«** 



c 



» (♦♦««*****«**»¥*«•* 



EXIT 

•4>«* *********** 
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SUIiHtlUTINE GNAOO 



LNTER 



* CALCULATE 

* NUM8FK OF 

* FtEHeNTS 



m***^* ********* 



************* 



EXIT 
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SlittHUOriMl CHSUi 



• • 

* CALCULA1E * 
« MOHSet) OF 4 

♦ UfcMtMS * 

* * 



* 4 

« S4J41TiIA4:T * 

« HAIRJCES * 

« ft>A-6 • 



EXIT 



c 



f^- 
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MVMMITIME CNMO 



•♦•••••••••—• 



! ao:^.: "K sis"^ 




•^ 



* • ••• 



• ta»***f set j»i 






• * • 

• iMm»«l* CLEMI OUTfUT 

• • • LOUTIM 



• mm 



JU, 



SCT !•! 



•iNCacasE K tv* 



■VMMM «1* 





I'M *»tl M* 

« • • 

• MO 



> m » n iii«iiii 

* • . mm* 
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SUeftOUTINE GNTftA 



ENTER 



* INITIALIZE * 

* POINTER AT 
vriRSr ELEMENT 

* OF OUTPUT 

* MATRIX * 



*«*•«•««***•••« 



HOVE ELEMENT* 



mm090*9* ******* 



* * *** 

♦NO * * 
LAST eLEHEMT *J* A2* 



( 



EXIT 
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SeURUbllhb oTv'RD 



CNTt.* 



• JhlllALIZt: * 

* PLlMtK AT 
♦FIHil tttHfcWT 

• Vf LtJlPUT 

♦ HATkJX ♦ 



««*«»*««*«*««*» 



• ♦ ♦ fKL« A AhC 4j* 

♦ ** ♦ « 



**♦♦•»«♦******* 



• KULTiPLV * 

* IkAWSPiii OF* 



«4S«4 «**«««*««« 



<**««««« 



L>.n 



»»«♦♦»♦**♦♦*♦*»*»*«♦**♦ 



4* **«««*««*«««* 

« ♦ ^** 

« iTfcP PClMtk* ♦ * 

* TC htXI *»♦♦)* Ai-* 
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SUeROUTlNt MAOO 



ENTER 



35 

••« • SET ARRAY « 

« « • l>01NT£K AT 

♦ B0****I*F1HST ELEHENT 

* • • FOR AUtH 
**« • MATRICES * 



* BOTH HATRlCeS •NO _ 

Of SAME ♦••»♦♦♦**•** 

* STORAGE *4QDE * 



»*»»♦* *• • ■ ■—**• 



••♦•*•••••••••• 



EXIT 



***«««««*«««««*«*****•' 



**4 • 

♦ * • 



«««*« «««••««*«« 



* • 

• calculate • 

• rTi»ji - • 
*Aii«Ji«eii,J)* 

* FOR ANV HSR * 



•NO 
LASf tLEMCNT •••♦•••• 



« • 

• CALCULATE • 

• HTCST - • 

• NSA*HSB ♦ 



EXIT 



«««««««««•«««« •«*««««« 



c 



• SET OUTPUT 

• SrORAfiE MODE 

• NSR-0 



• STEP PiilNTER* * * 

• TO NEXT •***1* Bl* 

• CLEMENT FOR ♦ ♦ • 
•BOTH MATRICES* *** 



*«••**• **«•• 



SET MSR - I 



•4«««***«*««**4*******' 



* * 35 

•NO • • 
NTEST GREATER •!* 80* 
THAN 2 ♦ • • 
•' « ••« 



«|« BO* 
« « « 
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SUBROUTINE NSUB 



ENTER 



••* * lUTAICES 



• ftOTH MATRICES •HO 

OF SAHE ••»••*••••• 

* STORAGE HOOE * 



*«*•*««*«*«*•«* 



MATRIX A 
«**«««** ******* 



** *« ****«4**««« 



• Bt*«*«|* 



** ♦ * * » ♦*» *** *** 



**II.j|-BiItJI* 

fM Mir («« • 



LAST ELEMENT »»»•«••»•••* 



*««•«*•***«•**««««**«•« 



«*•*! 



EXIT 



• * 

• SET OUTPUT 

• STORAGE HOOE 

• MSK-O 

• * 



• STEP POIMTEH* .• • 

• TO NEXT *M*I* SI* 

• ELEMENT FOR • • * 

*aoTH wkTmcES* *** 



•NO 
NTEST CKEATEII ••••♦♦»♦•♦••» 
riUN • 



• SET HSR - t 

« * 

*************** 



20 » » 

• • 35 

• * *** 

* *m * * 

* NTEST CREHTER *l* SO* 

• TH»N Z * * * 

* I * *** 

* ' * 
• * 

* YES 



•****«******«** 



• SET NSR - Z 



*i* 80* 

• * * 
*•* 
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^kUBROUTINE HPKO 



ENTER 



?0 

• • 

* 4 



* * 

PRODUCT * 
RUfKJ + ♦ 
*B(I*K)*A1J.II* 



i* 



•NO 

OI*GONAL BY •♦••*•**♦••♦♦ 

DIAGONAL • 



«*«*•*•«*«•**** 



* SPECIAL CASE* 

• FOfl DIAGONAL* 
«By OIACONAL ,* 
« R « A*B * 



****•*•**•••••• 



*•«*•««*«•***** 



• •«*«*««««*•«***•****•' 



•««•«**•******* 



•«****«*« «*•«** 



****•♦*•««** 



• • 

• STEP POINTER* 

* TO NEXT • 

♦ etewENT OT • 
*OUTPUT MATRIX* 



* •«* «* •« «*«« *«* 

* INITIALIZE * 

* POINTER AT ♦ 
•FIRST ElEHENT 

* OF OUTPUT * 

* MATRIX * 
mm************* 



***«««««•*•*«** 



•*««««•**•*•*** 



*************** 



•INCHEASE J BY****!* A7* 
« I * * * 



***«««***•*•*** 



♦*♦♦•♦»♦*•»•♦**•*♦♦♦♦•■ 



r 



•••••*•♦*♦♦•♦•• 



* A6***«l* SET J s 1 



•****« *««••*•*• 



**♦*♦••♦*•*•*•* 
••• * * 

* « « CLEAR OUTPUT 
« A7****l* ELEMENT 

* • 4 LOCATION 
*♦♦ • * 

••****«**•*•**« 



••«««*«**«***** 

« * 

* SET I == I 



*************** 
EXIT 

*************** 



********************** 

* 
* 



*************** 
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• INVERSION- * 

• -LOOP • 
••*«*********** 



•*••« *•••** **** 

* * 

* UPDATE • 
•I*IPIV-IPIV-MIN* 

*ANO IN0-INI>-1* 

* * 
****•«•*****•«* 



*•*♦**♦•••*•*** 

• * *•* 

• ADD 1 TO * • • 

• COUNTER K OF**** I C3* 

• RON-LOOP * * * 
« • **« 
•*******•«***•» 



******************* 



«*«*•«********* 



*************** 



********* ****** 

* * 

* ADD 1 TD * 

* COUNTER I OF* 

* WJLTIPLICA- • 

* TION-LOOP * 
**«««««••*«**** 



* SET * 

* A(1PIV)-0IN • 

* AND MIN-N • 

* * 
*************** 



*************** 



INVERSfON- 
• -I COP 

••***•**««•*«* 



IS I GREATER 

THAN N 

t * 



**«*****«*| 



I B6* 

• * 
**♦ 

• •• 

• • 

• B*** 



«•*••*••***•*** 

EXIT 
«*•***•*••••*** 
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!>JttKUUII\tr Slf^l. 



ENTert 



* t.LI'-il'Nfllt * 

• .»LL orH'.K • 



• SET KS TO 

♦ * 



LAST CULUHN 



•*•««•****«•**« 



•««****«•«****« 



* INITIALUE 

* COLUMN INOtX 

* * 



* a2****J*G£T PKECCKIho* 

• « • rttJU * 
• ♦* ♦ ♦ 



•♦• • • 

• ♦ * • 

• A3**«*l* CeT A COLUMN* 

• * * 1 



* SliLVt EUU 8V* 

♦ SUBSTITUTING* 

♦ PfttiVlUU&LV * 

* CUMPOT£0 • 

• VALUES ♦ 



•««**«4** 



* FIND MAXIMUM* 

* ELEMENT IN * 

* COLUMN • 

* * 

«•«««««*««****• 



* * 70 

« **» 

*NU * * 

ALL HUNS ♦»* rf^* 

PKOCESStU ♦ • * 



*««•*»«*« ***! 



• * 

« SET KS 10 I 

« 

« * 



EXIT 



••«•••*«««******«•*««* 



«*««**•******** 

EXIT 
•*•«*********•* 



*«*«*«•«*«***«* 



•••«•««•****«** 



SO V 

*************** b^ 

* * *** 
•OIVIOC BOW 8V* ♦ ♦ 

* LEADING •A**!* bU« 

• CQEFFICfENT • * ♦ 

• * *«* 
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SUBRGUriNES GtL& ANU DbtLb 



****•«***«***+« 



***«***«***«•** 



* * * 
AL •) A^* 

• * * 



IS COLUMN 
INTtkCHAhQt: 
NECESSAhV 



«•«•««*«« 



*•***«*** 



* • 

♦ INTtKCHANGC • 

* COLUMNS IN * 

♦ MATRIX A ♦ 

• * 



H *9******m*** ****** 



************ ••* 

* * 

* SET ERROA * 
|« PARAHEfEA * 

* 1£R— 1 * 

* • 



• ♦NG 

• IS TH£R£ LOSS * 

OF 

* SICNIflCANCE * 



* TKANSFUKH ♦ 
*PIVUT ROW AND* 

* INTFKCHANGt • 

* ROriS IN * 
« MATRIX A * 



*************** 



•«•••*••••**•*• 



*************** 

* * 

* set eKROR * 

* f>AftAM£TEK * 

* IER«K-1 • 

* * 
*************** 



*************** 

* * 

* SAVE CJLUHN ♦ 

* INTtrtCHANGb • 

* INFrjAHATION • 

* IN ACKtKI • 
*************** 



•«•«*•*•**•«•*««*** 



«( ***************** 



* SET ERROR * 

* PARAHETER • 

* ieR-0 ♦ 

* * 
*************** 



*************** 

* COMPUTE * 

* INT6RCHANGE * 

* INFORMATION * 

* FOR ROWS AND* 

* COLUMNS • 



*************** 

* * 

* TKANSF0K*1 • 

* RLST OF • 

* MATRIX A • 

* IF'JKMULA 61 * 



*************** 
*SEARCH MATRIX* 

* A FOR ♦ 

♦ ABSOLUTELY * 

* CREATEST « 

• ELEMENT * 



*************** 

* TRANSFORM • 
•PIWUT RQU AND* 

• INTERCHANGE * 

♦ ROWS IN • 

* MATRIX K * 



*************** 

• * 

• TRANSFORM ♦ 

• RtST OF ♦ 

• MATRIX R * 

• (FORMULA 71 • 
*************** 



• • 
«S£T PIV EQUAL* 

• TO AftSOLUTE • 
*VALU£ OF THIS* 

• ELEMENT * 
*************** 



♦NU ••♦ 
* • ♦ 
*l CO* 
* * 
*** 



*************** 

* SEARCH REST * 

* OF A FOR • 

* AbSiJLUTELy * 

* URtATtSr ♦ 

* ELEMENT ♦ 
*************** 



♦ COMPUTE * 
•TOLERANCE FOR* 

* TEST ON LOSS* 

♦ OF * 

* SIGNIFICANCE* 
*************** 



♦ »YtS 

• IS THERt ONLV 

ONE EQUATION 

* TO SOLVt 



*************** 
* * *** 

*StT PIV EOOAL* * * 

********* * TO AHSILUTE **♦*) A9* 

♦ *VALUE OF THIS* • ♦ 

* ♦ tLEMcK'T * •♦• 

• *************** 



*************** 

* INITIALISE * 

* COUNTER K-0 * 

* FOR * 

* ELIMINAnON « 

* LOtlP * 
*««*«**•****«*« 



L9 V 
*************** 

* BACK ♦ 

* SUBSTITOTIU.i* 

* AND bACK * 

* INT€l<CHANGt ♦ 

* rfO««UlAt 6}* 
*************** 



*t********^*0^0t*** 



*************** 



•A9**«I*AD0 I TU STEP****» BO* 

* * • COUNTER tC • « • 

♦♦ * ♦ ••* 

*************** 



i2 * 

V 
*************** 

txir 

******** 4****** 
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* • 

* PUT RESIDUES ♦ 
•INTO KIGHT HANO* 

• VECTOR * 

• • 



**f I******* 



««•******«• 






BIGHT HAND • 
VECTOR • 



APPLY ♦ 

PfHMjTATIONS Ti* 

«1SHT HBND • 

VECTOR * 



***«*02* ********* 

• S0LV6 THE * 

• SVSTe« USING ♦ 

• FACTaRED FORM * 

• • 



• IER*4 EXIT ♦ 



•*««*F1 ********** 



«•****«•*•«••««•« 



F2 *. 
,* EACH •. 
.♦ COMPONENT *. YES 
•. Of SOLUTION ,•-,.. 
«. IS .» 

• . «.* 

• • .* 
* NO 



* 

• I£R-0 EXIT 



ITtRATION *. 
STILL 
.CONVERGES.* 



.* I^ IT 
THc 2N0 
•.ITEPATIPN. 



• * 

• ItR«:* EXIT • 

• * 



HJ •- 

.*NCRH Of*. 
.•C3C?6CTICfN *. 
•• IS SHALL 
*. ENHUGH .* 



* lt.» = l EXIT 



•««*«HA*«* ******* 

* * 

* EPSI-NOSM Of * 
.x»C3RR£f II3N/NQ«H« 

* ri St^LUTIPN • 

* • 



» * 

* lt«=2 fKIT * 



fr. tir jnt.1l t\UlK,f 
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ll/Oi/67 



I'NTtR "AillR 



<■ ("QUI L 1 J'^AI I 



* 



• I ^ I 



««««<:«««*«*<: 



* SEARCH FOR 

* PItfJT IN l-T 

* COLUMN 



« caMPurE r-TH * 

♦COLUMN CF LCrf£«« 
* MdTRIX • 



• 1N0IC=3 



.«. IS DIAGONAL .* 
e. TEKM .* 



« « 



INfERCHANiE 

PIVOT ROW 4N0 

I-TH ROW 



COMPUTE I-Th 

kqw of uppek 

MJITRI X 
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SUBROUTINES HfCn AND DMFGA 



♦♦♦♦♦•••♦♦•*♦** 



ENTFR 



•* 



INITIALIZE « 

*R0***)« NCOL«H FOR » 

« « * ELI HI NATION * 

♦ ♦ * LOOP ♦ 



IS NUHBFR fJfVeS 
ROWS N UNO * 
NUHBtR nF • 

COLUHNS H • 

pnsinve • 



* *Nn 

IS 4es(Piv» * 

GReATEft THAN 

TOL • 



•«********C1* I 



•SET RESULTANT* 

* RANK IRANK • 

* COUAL TD -I * 

* • 



*^m** ********** 



**«••«*«**«**•« 



*************** 



m*^L* *********** 



*************** 

fXlT 
*************** 



*************** 

* * 

* INTERCHANGF • 

* ROWS Mi> * 

* COLUMNS IN • 

* MATRIX * * 
*************** 



******************* 



**t **************** 



*************** 



************** 



*************** 

* * 

* UPOATF ROW * 

* INntX VFCT.IR* 

* ANO COLIIHN ♦ 

* IND€X V-^CTlrt* 
*************** 



*************** 

* Ser Ut» THt^ * 

* MATRFX C ♦ 

* FXPRESSIN'; « 

* *",nvi * 

* OgcT^NTF^r, IPS* 
*************** 



*************** 



*************** 



IS IRAWt *♦♦***••••♦♦• 

♦ eOUAL TO N * 



*************** 
♦SEARCH MATRIX* 

• A FIX* ♦ 

• AeSOLUTFLY ♦ 

• GREATCST * 

• CLFMFNT ♦ 
*************** 



*************** 

* * 

* TRANSFORM ♦ 

* CURRENT * 

* SUBHATPfX • 
« * 
*************** 



*************** 

* SFT 0"* THF * 

* HATRIX H * 

* CONTAINING • 

* HON*ir.EN=0US * 

* SHLUTTHN ♦ 
*************** 



*^ ***************** 



*************** 

♦ * 
*SET PIV EQUAL* 

♦ in THIS ♦ 

♦ ELFHFNT * 

♦ * 

9************** 



*************** 

* * 

* SEARCH NFXT • 
•PIVOT ELEMENT* 
« PIV IN * 

* SURHATRiX * 



*************** 
EXIT 

*************** 



*************** 

* * 

* INITIAL rze ♦ 

* ROW INOEX * 

* VECTOR IRON * 

* * 
*************** 



*************** 



MDO N ro HCOi* 



*************** 



*************** 

* * 

* COHPUTF * 

* INTEKNM. •* 

* TOLERANCE ♦ 

* VAIUE TOL • 
*************** 



* «N0 *•• 
IS NCOL * • ♦ 

GREATER THAN *) Bl* 
N*H « « * 

• • ••• 
* • 

• * *** 

* YES • ♦ 

•**«***«**|* CI* 
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SUHKilUIINtS urLS A.HU UOLLi 



CNIt-t 



* • • * * « 

*rtO*) IS PIV tJUAL *l A2* 



• • * iAvr piv;;T * 

• * • MAT-^IX A ii^ ♦ 
♦♦ • VKCTOk AUX • 



•••*««*** 



• RdNSrORM • 

• PIV;! KO* A.^0* 

• iNTLKCHANGt * 
« AUUS IN • 

• MATRIX A • 



•• • • 

« « • SET CAAOH * 

**2«*«l* PARAMETER « 

• • • left — 1 ♦ 

•* * * 



• ♦NO 

• IS THfcBfc tUSS * 

Of 
« SIt;NIFICAIvCE * 



« • 

• SAVt COLUHN ♦ 

♦ INTfcRCHANGfc * 

* INFURHATIUN * 

• IN A<KfK) * 



tXII 



4 SET EKftOH « 

• PAftAMETEH * 

• ICR>K-1 * 

• • 



• TKANSFOKM • 

• RfcST OF • 
«1ATRIX A * 



« 
V 



•( ««*****4**««**««* 



S£T ttMon. 

PARAHCTER 

IER-0 



• COMPUTE ♦ 
« IMTERCHANGk * 

• INFONHATIUN « 

• « 



• • 
« TKANSfORM * 

• HtST OF • 

• -iATAIX R ♦ 

• (FOkMULA 71 • 



• OETERNINE • 

• ABSOLUTELY • 
« GREATEST • 

• DIAGONAL • 

• ELEMENT IN A* 



4 TRANSfOKH « 
•PIVOT MU ANO* 

• INTERCHANGE * 

• ROUS IN * 

• MATRIX « • 



•SEARCH ABSOL.* 

• GREATEST ♦ 
« 01 AGONAL • 

• ELtMtHT IN • 

• REST OF A • 



• * 

•SET Piy EOUAL* 



•NO ••* 

• • * 
• J CO* 

• • 



« « •*• 

•SET PIV fcUUAL* • • 

• TU AbSOLUTE •••♦J A9» 

• VALUE OF 1 His* • • 

• ELEMENT • •*• 
•«•«*••«••••••• 



• COMPUTE • 
•TOLEftANCE FOR* 

• TEST ON LOSS^ 

• OF • 

• SIGNIFICANCE* 



♦ •YES 

IS THERE ONLY • 



•••«••••««•* •«* 

• INITIALIZE • 

• COUNTER K-C • 

• FOR * 

• ELIMINATION • 

• LOOP • 



JO V 

•••••««•••••• 

• BACO. * 

• SUBSTITUTION^ 

• AND MACK • 

• INTERCHANGE « 

• (FORMULAE: SM 



•(«••••••«*«•••••«• 



• A««M|*AOO 1 TO STEP*«*«I BO^ 

• • • COUNTER K * • • 
•• • • *•• 



EXIT 
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<^UilRUUriN(S i;tL3 «NU UlitLO 



ENTER 



•• ♦ CUNPUIF • 

* « •TULERANCe FU<4* 
>UO»«l* I£ST UN lOSi* 

• • • Of • 
** • SlfiNIFIMNbE* 



•• • KCUULE 1*1 KUT* 

« * « RUM AI40 • 

*CU***I* imtRCMANCt • 

* • * RIMS IN • 

»♦ • H«TKIX A • 



ME INPUT 

PARAn£TEftS 

UMRECT 



• SET ERROft * 
I* PMAHETER • 

• lER— 1 • 

• * 



EXIT 



•*••«««*«*•««••••*< 



* JkllliuuE • 

•OECOMPCSITIUN* 

* LOOP «N0 SET* 

* COUNTER K>0 * 



••••••««• 



•«««««*•••«•«*« 



UPDATE 
ADDRESS 
PAftANETERS 



♦ •ves 

I $ K-TH RUW • 
LAST Raw TU 

•TRANSFOftHEO* 



••••*•••• 



•TRANSrORM THE* 
•OTHER ROUS IN* 

* NATHICES A • 

• AND R * 
•«*«**•*••«•*•• 



«(•***«*«•******•** 



• UPDATE • 

• ADDRESS * 

• PARANETERS * 
« * 
*************** 



*••*•*••*•*••** 

• oescrTSinc * 
« natrik a and* 

»*•***•*••*•*•• 



«•••••••««**•*• 



COIUIIM OF 

• NATRK A • 

•**•••*••*••«*« 



HAS NATR 
LONER 
COOIACONALS 



14 

•NO ••• 
A * • * 

•I ao* 



• * 47 


• * 


« • 


« • 


•YES •♦* 


• tVES 




* MS INPUT * 


IS PIVOT •! A2* 


• MATRIX A OF 


EQUAL TO * * * 


* DIAGONAL FORH • 



********* 



*************** 

*MOVC ELFMfNIS* 

• OF MATRIX A • 

* 8ACKHAH0 * 



0£rCRMIME * 
AaSOLUTELV « 

WEATESr • 

6LCMENT IM • 
MAIKIX A « 



• *Na 

« IS THEItE LOSS • 

CF ***•**' 

* SICMIFICAHCIf * 



• Y€S 

« 
* 



*************** 

* • 

• Sti CRrttlH • 
« MRAftCTEft • 

• IER-K-1 • 

* * 
*************** 



* * 
*dACKSlMSflTU-« 

* TE ELENCNr^ * 

* Of HATRU ft • 



*( ***************** 



Ab • 

V 

tXIT 

*************** 



INSgRT 

KIAMCULA 



26 V 

*****9********9 

• NEOUCE P1</UT« *«* 

• KOU AKO • • • 

• INTEHCHANlie *«««| CO* 

• KOi^S IM • • « 

• NAritlX H • *** 
*************** 



in 



SUeXOUTIMES HTOS AND OHTOS 



ENTCB ♦*«*••♦*♦•«•♦• I 

* 



• tNC 

AftE * 

01 HENS IONS •#«*•*»♦♦••♦••♦•• 

H,N POSITIVE • 



>«•**** 



m ****** ******** 

• INITIALISE ♦ 

• DIVISION • 

• PROCESS FROM* 

• LEFT HAND • 

• SIDE * 
*************** 



***************** 



*^»*m************************* 



« * *************** 

•NO * INITIALIZE * 

* • DIVISION • 

POSITIVE ^ ♦ I TPAH^SEO • 

^ « *************** 



V 



ABSOLUTELY 

GfteATER THAN 

* 3 * 



* • SET ERItQft * 

♦ ♦♦*♦**•••••••#♦•*•♦•♦•*••*••••••••**•**•**••*•♦ I* PARAMETER • 

* EQUAL TO -I * 

« * 

*************** 



.Y|S 



• *«•««« 44 *ft5*< 



• INITIALIZE 4 
•14 SOLUTION Of 4 

• T»X.« S€T ICR* 
4 EOUAL 1 4 

444 44 44 4*44 4444 



444444444444444 

EXIT 
444444444444444 



.44444444444444 444444444444444 

* INITIHIK » J J 

iJiriEi ««ull4 4 DIVISION 4 



444444444444444 



444 444444444444 



4444444444444444444444444444444444 

4 

V . .. • 



IS TRIANGULAR 

FACTOR 

SINCULAR 



************* i 



*************** 

* 
* 

***************i* 
* 
* 
*************** 



SET ERROR 
PARAMETER 
EQUAL TO 



• SET ERROR • 

• PARAMETER * 

* EQUAL TO I • 

* * 
*************** 



IS lOP NOT 
* EQUAL -3 OR 3 



c *************** 

^ * ♦ 

*************%* EXIT * 



*************** 



* •** 

• NO • • 



*************** 
EXIT 

*************** 



SUeROUTINES HLSS tNO OHISS 



• ♦••••♦*•**•••* ♦ •NO i * Z I 

♦ » ♦ • • • • 
ENTER ........«....! IS DIUCNSION ....«....•••♦...•••♦••♦ J''".**'l'itilikBlSl Jit t 

* • N POSITIVE • • • • 1 IMNSPWI'XXJ 

•♦•• ' . . * ; ;*«••*••••***«* 



«<Sk IMNK »»«»«..»...•«.«.»•*». |. • IS ISJ ESUAL •! M • 

POSITIVE . : • . ' , * ;..5 



V « • 24 

res • * 

* * SFT FRROII 4 , . ~ 

,n,«.*>*>»«**|« PMMETEli • »02^« ^ IS ISjf BIU«l •! »l 



•TES •••• 
* • • 

eou*i TO -I . ;. • • . J . • • 2» 



• SET e»RO« * »«»».»•»••••»*• 

• PWMETEIt T0» r ,_„ 

• l€BO • « MIT 

• INITIALIZE • • „.„.„.„..„ 

• PEHWTATISn • ».«..««»»»♦»»*« 



• 


« 


* 

* 




16 






«VE$ 


*•«« 






• 


• 


• 


s 


ISW 


COUAL 


•1 


A2 • 




1 


• 


• 


2* 



* HANOSIOE * 

* • 



pMwT«flOH * ....««.....|J EXIT 

PERFORMED ^ • • „„„,»„„,„ 



• SET SWITCH * 

* ISW EOUAL TO* 

: ' : 



*«•*•«•««**«««*•«**•*«•* 



IS DEFECT *••♦*••♦•♦•*••♦•♦••♦*«)♦ 

MUAL KfO * * 



• ♦ 

* • « 

• • 28 V IT 

•♦* ♦ * • •¥£$ • INITULIZE • ♦♦*♦ 

• « * * • * ♦ DIVISION Of • • * 
*A9 ••>♦ CALCULATE ••••*•*•♦♦••••*• I IS ISW EQUAL ♦♦♦♦♦♦♦♦••♦••••|*J(l SET SWITCH****! A3 * 

• * • Kl*)Cl-fU*X2 » • 2 • • ISW EQUAL TO* • ?• 
*«* * * * * • 2 • **** 

* « *** 
• HQ * ♦ 

* * 
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SUBDOUTINES MLSS AND DNISS 



•«•••**««***«•• 



«*0 "I* Of iLeKfUTS ••••I 02 • 
• • • Of « • J^l* 



••• * • 

•*i **l* ISM caiWL TO* 

• * • 9 • 

••• • • 



••••«••** 



• • * xT 




* SINCULAR • 






* E«UAl TO 1 * 



r 

V 



• IS IS* fOIHU. ^ 'J 00^» 



.V|S-*., 

• I M * 



. ^ U II| MUM. , .J^ 
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SUSmuTINCS MCHe MO CMCMe 



CNTW 



FtCTMtltTTOM 



««••••«•«« 



LESS THAN 

Jew • 



«•*««*« 4#*«******«*»«*| 



•«*•«***# 



* SET flC-«UO*l* 

* » 

* • 






IS giMENSION • * • * SET 

M LESS THAN *M******n4 *•)• hui 




««***X****« 



• INITIALIIE • 
*ll(W-iNiEX K-1* 



•«*****Mi 



0* 

-I* 

lEK— 1 



•«M***M«***C4 )• EXIT 

• * * • 

• •** MM*********** 
••••••«*««««*** 




**<****M ****** 

* tI" • 

****••*«***«*** 



• SET vr • 

* moH-NWi AND* 

• ERROK * 

* rAKAHETEII * 
» lER-O • 
«****«**M***** 




THIS IS 

WANCH TO 

FACT (MfZAT ION 

OF MATAIX A 




• YES ♦ * 

*•********! « Cl • 
♦ 2* 

2A 



CIVISION OF A 
TRANSroSEIAI 



* IS lENtTH OF » • * 



•MMM*«*M*M 

: Fi8SPi.s: 
LgsSt!*4 



************ *M 

* *« 

LI2E * * 

- |nte«****i 

. OA JNM* * * 

FAOOOCT l-l * •• 

*************** 



IMITIAUIE 
STEF-CoiNTE 

FOA JMK* 
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SUemUTINES mchs amo ohchb 



I* • • 

is them lossno 

oTgits ih • • 

•rao ic and suh • 



STAUT 

DIVISION OF 

HATRIX « BY 

TRANSPOSE OF 

MATRIX A 



- - •tl-tl-tl6.AM0« 





•«•* • •ves •*** 


. 


** » * IS ABSOLUTE • » » 
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• INITIALIZE • 
» STEP-COUMTi** 

• FOR INNER • 



ISTRACT 1 • 



' RON' 



•<#IMTE LENGTH* 

• OF INNER * 

• PRODUCT LNAX* 

• • 



•*•«•««•*•••••• 



• SET UP • 

*ROU-INDEX K-H* 



*************** 
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DIVISION ay 

TRANSPOSE OF 
MATRIX A 



«•*•••***«**••• 
• • • • 

• • • INITIALIIE • 
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•«*«*•««« 




.r|S 

IS K S90M. TO «»«... 



ESUM. 
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I* • 10 * * • 



10 
* TES 
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SUBROUTINES PRBH AND OPRBH 



♦ • 

• ♦Str NUMBER OF* 
•«1* CALCULATED • 

• *A00rS lR-lR-1* 

* * 



«• * « 

« • *SET NUMBER OF* 

*BO«**l« CAICULATEO • 
♦ • •ROOTS IR-IR-J* 



* « * SHIFT BACK * 
♦Al ♦•!• £L6HENTS OF • 
« * *VECTOR POL BV* 
•** * I * 



****••****••••• 



• **«*****«**• 



*****•*•••*•*•« 



•«««••••***«*«* 



•«« • * 

• ♦ * SET ERROR • 
•02 **1* RARANETER • 

* • • IER-2 • 

*** • • 

* •• «* *• *••« •• •• 



«•*••••••••«••* 



POL 

• • 

«•«»«*•• ••• *** • 



• • 
♦SET NUMBER OF* 

• CALCULATED ♦ 

• ROOTS IR>0 • 

• * 



• STORE AflSOL.* 
•CREATEST REL.* 
*OIfFERENCE IN* 

• RR<IR*lt * 
•«• *« *4 4 4«« •** • 



«***««**«**•«•* 



• •• •• «««•*« •*** 



V 



IS lER EOUU •••»••«« 
TO zeRo • 



« *N0 

• roOK JtCCUItCY * 

Of RESULTING »»•»..•. 

• COEFFICIENTS * 



• SET ERROR * 

• PARtHETER * 

• lER — 1 • 

• « 
♦«•♦♦♦••••*•*•♦ 



*(••*****•**«*****< 



41 • 

Htf ********* 

EXIT 
*************** 



• * 

• SET ERROR • 
I* PARAMETER •• 

• IER»3 • 

• • 
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SUBfiauTINES POFS MO OPQFB 



ENTEK 



•• • OEflNE • 

• * • TOIEMNCES • 
•B0***|* EPS AND E«l« 

♦ ♦ • FOP TEST OF • 
• • * ACCUKIIGV • 



• • • COHPUTE NOftH* 
•CO»»»l»M OF WOIFIEO* 

• • ♦ LINEAR • 
•* • KENtlNOER • 



28 
•• • SET HH E0O«t« 

• * * NAXinUH OF * 
tO-^tttf «BSfllUTE • 

• • * VALUES OF • 
•• • A1.B1.C. • 



••****«•««««»«• 



*««*«**• ««••*** 



• INITIALIZE » 

• ITERATION • 

• COUNTER L-0 * 
AND SNITCH • 

ll-S • 



It-i 

• *«*# 



•NO ••♦• 



IS SNITCH LL 
• EQUAL TO ZERO * 



*l 6*^ * 
• 7* 



• OETEimiNE * 
•ACTUAL NUMBER* 

* OF * 
•COEFFICIENTS.* 



• SET OP ♦ 

• PARAHCTERS • 

• Be * 

• •••••«4******* 



19 



• ADO 1 TO • 
» ITERATION • 

• COUNTER I • 

• * 
••••••••••••*« * 



• NORMALIZE • 

• COEFFItlFNTS* 

• Al.BliCI.A.a* 
•DIVIOIM. MITH* 

•••«••«••••*••• 



*l A6 
• 2' 



•COEFFICIENTSYES •••• 
OF LINEAR * • • 
REMAINDER *l A4 • 

SUFFlCIENTtV • • 2* 

♦ SMALL • •»•• 



•****•*•*•••*•* 

* « 

• CALCULATE • 

* DETERMINANT • 
•H>Al*CI-ai*ei« 

• • 
**•*••*•♦••••♦♦ 



• DIVIDE • 
« COEFFICIENTS* 

* »T LAST • 
•NONZEilO TERM.* 

I.E. CIJI * 



*«****«*•**** ** 

* INITIALIZE * 

• HININUH NORM* 

• cc • 

* • 
••*•* «««••**•«* 



• 'YES **•• 


• PREPARE TEST* 




• OF NORH OF • 


OVERFLOH IN *l AB * 


* MODIFIED * 


NORMALIZATION * * Z* 


• LINEAR * 


• * ♦♦•* 


* REMAINDER * 




*************** 



* 'YES 

IS NORM H * 
CREATES THAN 
MINIMUM NORM 



22 V 

**«****««***•** 

• UPDATE * 

• HINTHUM NORM* 

• CC AND • 

• PARAMETERS ♦ 
•AA.BS ,00 1,002* 
•*****•**••**•* 



•YES •*•* 



• •••***«< 



*( *******••«••**** 



DETERMINANT H 
* EQUAL TO ZERO • 



30 V 

*•*******«****« 

* COMPUTE * 

* REFINEMENTS • 

* OF 01 AND 02* 

* AND REFINED • 

* 01. 02 * 
*************** 



41 

*YBS ***• 
* • • 

*1 AB 

LINEAR * * A' 

* **«* 

* 
* 
* NO 






l«****** 



23 • ♦ 29 

• * 

IS COUNTER L * * « 

WEATpt THAN •» 0^* 

Ll« * • * 

* « *•* 



* *NG *««• 

IS THERE * • * 
5ATISF«T0«V •) AO • 
ACCURACY * ♦ 2« 



IS POLYNOMIAL *l BO* 
QUADRATIC * * * 
• * «*« 



444*4 •#«*•« *•*« 



HAS NORN H 
REMARKABLY 
DECREASED 



44444*444444*4* 



4CHANGE StiITCH4*44) B6* 
4 Lt TO 1 4 4 * 

4 • 4«* 

444444444444444 



30 V 

44***4****4*44* 

♦ SET RESULTS * 

• 0I3I-O, * 
» OI«l"0 * 

***•**•*•«•**** 



* DIVIDE * 

* POLYNOMIAL • 

* TKICE BY * 

* CUAOBATIC • 

* DIVISOR • 
***** ********** 



25 * • « 

• « 

* *YES **•* 
MAS X*)( • * • 

INITIAL eUESS •! t" * 

OF ouioRiric * * z* 

• FACTOR • **•• 



*************** 

EXIT 
*************** 



*l 80 * 
• Z* 

**** 



EUESS OF 

QUADRATIC 

FACTOR 



•*********| * CO* 
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sutaouTiNts rva mo DPSfs 



42 



*M * 



IS COWNTER L . 'I •» • 



•to *•>* UiiMm* * 

««• • 



*i 



<•*«**•««••' 






• ves 



•***«*««*MMM 



• ^uiu wwj • •*»»>•*♦♦••♦ 

• * ■-" * 

• YES 
••*«4**« 



■•* 
12 



•<♦♦♦♦••••••••• 



••* • SEI «{ 



*♦♦•♦•••♦*♦•••♦ 



• SET EWUMt • 



EXIT 



*:i* 



* * « SET EMOIt « 
•Aa *»»« rARANETfK • 

* • * IEIt«-2 * 

•«* « • 

• ** «• ********•« 



EXIT 
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SUUOUTINES CNP AND OCHP 



• SCT Till 

* BHMl. TO I 



IS moeR 

OtClTE* TM 






*************** 

* 

»• EX IT 

* 

*************** 



EiiWTIOM • 



* IS I SUfATEK • * • 

• atoti * * * 



EXIT 



151 



SMHOUTIIIES CNPS «M0 OCNPS 



{NTCD 



OlAfNSliN ' •••»•••«••«*• l« EXIT 



•J ™L 

-IMCNSION 
POSITIVE 



* • 

• I * * * 

* * ! 



«••••«• 



♦ SET «II««X*I ♦ • EXIT 



• INITIAillC « 
•HO»Q Nl-0 SET* 

• I'M • 

• • 



« • •HCpLACE Mi tV* 
4«9«**l« HI AND Hi BY* 



• COtMTEP I (T» 

• I • 



• •Y£S »•• 

• * • • 
IS ■ SKfMW •! *5J 

• TMW • • „• 



• • 

• set v«o-5« • 

• • 



(III 
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SUBOOUTINES TCNf iMO OTCMT 



ENTER 






IS THE 

OIIIEMSKM 

raSITIWE 



« RCfLACC * 

* coeeficRnt-* 

• VfCTOK OF • 



« EQUAL cai • 



• coeBicilHi-J 

• VECTm OF « 

• TfJ-i.XI BV • 



« SCT rOLdI • 






i cfcsfil's* 



• SET FO^I?! ' 






A4 4 4iAi#v4## 4i#A 



IS rot 

IMEHSION 

tOUM. 2 



* 
•I* EXIT 



IS J EOUAL TO (I CO* 
H • * • 

• • *** 



• yETM F0« 

• iTS.Ii AND 



V 

e«iT 



• INITIAL IZE • 

• ITESATION ** 

• COUNTEK J-2 • 

• • 
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SUBROUTINCS CSP MO OCSP 



ENTER 



»mm—*m*******m 



•HO *«««*««*«4«**«« 

***4«**««*«**| • EX IT 

4 



• SET ¥I2J • 

• EOUAl TO * 
« 2*Alt6UftENT-|* 

• * 



•••**«***•«**)• EXIT 

««««*«««*««**4« 



( 



•«•«••••••*•••• 



««««*««**4««*** 



«M***|* SI€P INDEX I* 



*•*••*•**•*••** 



* CMCU4 

• IIECURR 



■ FftOM 

KRENCf 

EOIMTION 



• (NO **• 

IS I atstta » • • 

THiN CIVE« •■ U* 

anoH tit 



EXIT 
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SUBDOUTINES CSPS WD OCSrS 



*•«« ******♦••*« 



••«**•*•«•** •** 



POSITIVE * * 



««*«**♦*•*•«««• 



•*•**•*****•*** 



f , ♦♦♦*••♦♦•*♦•••• 

• *"*' ! * 

DIMNsVoN ' ****»**»***——>* SET JCSULI ♦ 
GREHTeR rioj • J T-C<ll ♦ 



EXIT 



• AND THEH • * f«" 

• ARC-«tC*WC • 

• • 



* * 

* INITIALIZE • 
•H0>0 M1"0 SCT» 

• l-N * 

• * 



•* • • 

« « WEPLACE H2 BY* 
*A9«*«)« HI AND HI BT* 

• « * HO • 



•••♦**•«««*•«*« 

• SET WAIIG* ♦ 
« H1-H2»CII> * 

• DECREASE • 

• COUNTER I BY* 
» 1 ♦ 



I GReATCR •! AS* 



««««*«*•••««••• 



EXIT 
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SUSAOUTINES rCSP AND OTCSP 



ENrEA 



•C0*««1« CO€FFICIENr-« 

• * • VECTOR OF • 

•* • TS<J,X) • 



IS rHE 

OIHeNSION 

POSITIVF 



•NO •*««•««««««*««« 



• REPLACE • 

• COeFFIClENT-* 
VECTOR 



•««••*••••««•«« 



IS THE 

OIHENSION 

CaUAL I 



••*«•••••••«*•• I « 



£OUU. ciil • 



VECruK 

• TSU-l,Xt BY* 

• TS(J.XI • 



* • 

• SET raj.<ii 



1} * 
>! * 



••••«*««••*«••• 



CXIT 



• • 

• STEP • 

• ITERATION * 

• COUNTER J e«* 

• I » 
*************** 



*************** 

* SET I>at.l2l • 

* ESUU. • 

* CI2I«2» • 



: ThH-" : 

* TSIJ-1(I(I TO* 

* COEFFICIENT-* 

* VECTOR POl • 

*************** 



V 



IS THE 

DIMENSION 
EQUAL 2 



*•«•••«**•«*•)• EXIT 



• « « 

IS J EQUAL TO «) CO* 
N ♦ • • 
• • ••« 



• V£S 



* INITIALIZE • 

* COEFFICIENT-* 

* VEC TOR • 
(O.XI AND • 

TSh.XI * 



••••*«••*•••««« 



*«••*«•««••*••« 



• * 

• INITIALIZE * 

• ITERATION ** 

• COUNTER J-2 • 
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SUSMMniNCS HCP «N0 OHfr 



MT VttI 
MUM. to 1 




I* EXIT 



.T|, 

IS THE DHaEII 
* MUM. TO I * 



: m'Aii't 




'fni.'SS? 


« 


• 
•1 


* 


IWOCIl 


* 


• 


* 






•«* 


• * 








* « 








• res 
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SUBHOUIINES HEPS *N0 OHErS 



CNTEK 



r 



IS THE 

OHKNJIOHI 

POSITIVE 



* 
•*«•*«♦«•«*** I* EXIT 

* 



It THE • 

OIHfNSIOM 



******•*«««•••« 



* INITIM. 

• MO"-' - 
•Hl-MI 



IM.IIE • 
IIOt>l • 
.«• MT* 



•*9**aM SET 

•_• • l«-X»H»- 



MEPLitCE NO (T* 

• MI, Ml (Y * 

7.^!ff.tT?.!i7 



MT 
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SUSMIUTINES THEP UNO OTHEP 



ENTER 



• » • CM.CUL«rE • 
•CO*** I* COEf*ICI£HT-» 

• * * VECTOR OF * 
HIJ.XI * 

- ' - »i*****m 



■ S THE 

DIMENSION 

POSITIVE 



»l« EXIT 

•*••«•••*••*•*• 



* RErLACE • 

• COEWICIENT-* 

* VECTOR OF * 
» Hij-2>xi ir * 

• Hij-i,«i • 

********* 



IS THE 

dihension 



******«•«•«««•* 



•******«4*«**** 



SET POtlll 
EQUAL cut 



SET PO 



EXIT 



• REPLACE • 

• COEFFICIENT-* 

• VECTOR OF • 

r«*«Z****««* 



* STEP » 

* ITERATION * 

* COUNTE* J BT* 



** **»♦»*—* ♦*• 

• * 

* SET ran2t * 

* EOUAL • 

• CI2I*2*A • 



ADD CIJI • 

TINES * 

HtJ-l.XI TO * 

CaEFFICIEHT-* 

»»*5*** ******* 



IS THE 

OIHENSION 

EQUAL 2 



* 
*l* EXIT 

*•••*••*««««*** 



• vecTO* FOR • 

• HlOtXI AND • 

• HU.XI • 
«*•«*••*•••«•*• 



•***•*****• 
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suwunniHES lap mo olat 



^SL^Wt 



• *N0 

IS THE OKOHI • 
OtCATEX TIMM 



iKI^T 



( 



• JNiriAtIZE • 

• INDEX l«2 • 




t two • 

IS I SKCAIEII * • 
THilN CIVEN *l 

oriocx • * 



MIT 
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SUUOUTIMS LAPS iMO OCIK 



•••••♦♦**»i I ii» 



•««««*«««*4««** 



IS THE 
OIMCNSION 

rasiiivc 



**•*******«*••* 



••**••**••*••*• 



IS THE • 

DINEHSION 

GtCArea imm • 
* I • 



•I* EXIl 



• * 

• INITIIU.IZE * 

• MO-L(OvXl * 
•Ml-ill.Xl SET* 
•1-2 «N0 r-i-x 



** • • 

• • • SET • 

» • * IT*Hl-Hal/i • 

*• • • 



•««•«**«««•««•« 



*«*««*«***44*«* 



• too Ctll'HO • 
•10 » AND STEP* 

• COUNTER I BV* 

• 1 • 



•NO ••• 
• « • 
• » «5» 
• • 
••• 



•••••••*•••••«• 

EXIT 
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SUBROUTINES TLAP AND OTLAP 



ENTER 



* • 

* CALCUtATE 
COEFFICI" " 

VECTOR _. 

* L(J,X) * 



IS THE 

OinENSION 

POSITIVE 



***•««*«*•««# I • 



EXIT 



* REPLACE ♦ 

* coeffIcUnt-* 

* VECTOR OF • 

* LlJ-2«xi ay ♦ 

* L1J-1.X) • 



««*•«»••••««*«* 



IS IMF 

Df HENSION 
BQISAL 1 



**«*«*44**«*«** 



* REPLACE * 

♦ COEFFICIENT-* 

• VECTOR OF * 

• LU-IiX) 8Y « 

* LIJ.Xl • 



• • 

• SET POHl> * 

• FQUAL • 

• C(lt«CI2» * 

• -B*CI21 • 



exiT 



* • 

* STEP * 

* ITERATION * 

* COUNTER J BY* 

* 1 • 



«««««••*«****** 



**«•«**«•***••• 



• ADD CIJ) • 

• TIMES • 

• LlJ-l.Xl TO • 

• COEFFICIENT-* 

• VECTOR POL • 
*************** 



c 



*YES 
IS IME • 

DIMENSION »««»•«*•••«» 

EOIML 2 * 



EXIT 



• WO ««* 

IS J Eau«L TO •! CO* 

N ... 

. . *** 



• YES 



..«.♦*♦...♦.•.. 

• (NITIAIIZE • 

• INEFFICIENT • 

• VECTOR FOR • 
. LIO.XI AND • 

• lIl.Xl » 



EXIT 



. . ... 

• INITIALIZE » . ♦ 
. ITERATION •»•«) C0« 
» COUNTER J«2 » • » 

* * ... 
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SUBROUTINES LEP ANO OLEP 



•««*••**«••««*« 



COUiU. 



iCAIEK THM) 




EXIT 




IS THE OROER 

eawu. TO 1 ' 



I E»ll 



•««««««•«•«••«• 



*M*«*I* STEP INDEX I* 



CALCIM.1TE • 
riii niaH • 

li5*TfgS' : 
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SUMOUTINES lers «HD IN.EPS 



«*•«««««••••«•• 



IS iHe 



**«**«***««««|« EXIT 



plHEMSION 

oie 



•«•*••••••••* I • 



*••••«**«•••••• 



1-PI1<X| SE 



• * • ttt-MJ-HOMU-* 

•• • 7h - 



*************** 

* • 

•KCTLME MO n* 

* HI «N0 HI ir* 

* H2 • 






IS J SKEATEII •{ AS* 
IMUI N • * • 
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SUMOUTINES ILtr MO OTlfr 



•C0***|* COE^Y^IENT-* 

• * ♦ VECTOR OF * 

•* • rtJ.XI • 



• IS THE • • • 



* coeM[cRnt-« 

• VCCTIM Of • 



Diil«tiaN 
• twSi. 1 • 




• set roLUi • 



• ca^VcflNT-* 



• • 

• ST» • 

• ITEMflON • 

• COUHTEll J «r* 

• 1 • 



^SLr&iiii 



•««««*< 



• PU-t,>l TO • 

• VECTOR fOL • 



• •YCS 

** DlilN»&. * 

• COUM. 2 • 



IS J EOIMI. TO •■ CO* 
II • • « 

• 4 «*• 



iiii lnnn 



« 
♦ um i l i 




af • • • 

J!!, r-i "t 
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suaRouTiNES mm urn okthi 



••«««*«««««««•• 



•*««••«• *«««««4 



* * * COMPUTE NEW • 
*aO*«*l* KEFINENENI A* 

• • * AHO ITEKATEO* 
... ^ 



**•«***#*«*«*** 



SET ERKOK 

PMAHEtER 

IER>0 



•••«*«•«••**•** 



*************** 



*+«*•« ••*••«*•« 






• INITIALIZE • 

• REFINEMENT * 

• «<X->ST INO * 

• SET •— * • 



***********mm** 



*************** 



• INITIALIZE • 

• V*L'X-FCT(«i* 

• «NO • 

• STEI>-CauNTER« 
l»0 



*• * * 

• * • *00 1 TO * 
•«6***l* STEP-COUNTER* 

• * * I • 



• • 

• SET ERROR « 

• PARAMETER * 

• IER-1 * 

• • 



•{•**««*«*«***«**•< 



•••*■• EXIT 



« COMPUTE * 
* OCNONINATQR • 
•OF REFINEMENT* 



IS 

DENOMINATOR 
EQUAL TO 



•••*•*••*•»•*•*>« 



SET ERROR 

PARAMETER 

IER-2 



*«*•*****••*••• 
EXIT 



• **« 

» YES • • 

*********«|* BO* 
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SUMOUTINES mill MO ORTNI 










,imm, 



* FUNCflai • 



• • • «D0 I TO • 
*a2 Ml* STEr-CaUMTE«* 

« • • I • 



* INITltLIK • 

* stcr-COUNTEIt* 

• MOP... • 



• * • INTE)tCH«NCC * 
*COM*|* LEFT MO • 

* * * ItlSHT BOUNOS* 



♦• • •MO ••* 

* * * IHOULO • * * 
•Cl*l SISCCTION BE *l DO* 

• • • REPEATED • • • 



• * •coMfUTE OX sr* 

•00***l* INVERSE • 

• • • PiWAWtlC * 
•• •INTERPgLATION* 

•*******«3***«* 



*» ** » **»» » ♦ * ♦ <♦ 



• YES 



* uroRTE n-« • 

* AMD FR-r FOR* 
•MEN 81 SECT ION* 

* * 



vtolerSnci 






* CEMeRATE • 

* R€FINEO X AS* 

• X-XL-OX * 

• • 



^m\h 



• i • • 
IS KT<St.l-0 •• C«* 

* • i • 



• MITIRLIU • 



•. ^miiit *• 



• STEf^COUMTER* 

- K • 



'■*«****■»* 

• * 

• GOMPUTC m * 

• X>iJlC*JIR}/Z • 



>•«****«• 



SRrlSFACTORY 

I ACCURACY SV ' 

* BISECTION * 



• (NO 

IS K fl*'*'- "> 



...JU.*..*. 



1 A4« 



* S£T EUUli 



* * 1* 

• *Y|S ••* 

IS FCTIXI-D •» C»* 

• * * 

• * *** 



*u 

* LOOR • 

- - - -«•***•*••• ' 



*YES **• 



• SAtlSFACTORV *i 09* 

* ACCURACY BY • * * 
INTERFOLATiON •»• 



• •YES «•» 

• • • 
IS FCT<«RI-0 •( C«* 



*YES *•* 
IS FCT(XI-0 *l C»* 




• CNO ••• 

• SHOULD XL AMD • * • 

XR BE »t CI* 

• INTERCIMNCCO • • * 



** ISFR '^ 

ABSOLUTELY 

• SREATER THAN • 

* PL » 



* * 

««••*•««««•**«« 



*( **•««*«**•«*««*•< 



,Y|S «**, 



23 



* • ••• 
*SeT NEM RISHT* • * 

• BOUND XR-X *>**) B2« 

♦ mn-f * * ♦ 

• * •** 

*•*«*••«•««***• 



• •YfS 

* tt*f FL ANO * 
^ (R THE SAl3 



•I «K^« f 



»•»• EXIT 

- « 
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SUMOUTINES RTHI AND ORTMI 



• • • UT HCN UrT* 



i^iisrr-iy 



V 
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SIMKOUIIHES RINI «N0 DRINI 



•bo***i*iiefInewnt ox* 
•• • mcunent x • 




• ccmruTE • 
FUNCTION FIXI 
AMO 
0|Hj(fT!VE . 



«***•*««•«««••• 








• »rES 

■S THEAE • 

SATISFACIORT * 

ACCUItACT • 



* 1-0 



*«a ••• 
• • • 

IQ »l AS* 

« * * 

• •** 



1 TO 

cauNii 



•* 




•J 




SCI ERROR 






PARAMETER 






IER-1 








•• 



•I *«««««**«»******l 



••••*•*««•«***• 

EXIT 
•***«•**«•*«*** 



IS ITERATION •! BO* 
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SMIIOUIINES FNFr MO OFNFF 



•*•«««••••«*•« • 



ROOUCTt 
RtX - 




• I *9 • 



*i 



IS OIRECTION •) *0 • 
VECTOK Smu. * * 2« 



COtMnH Tl 



• •NO 

IS HMUM/CNIM • 
MEnER TIMN 
EPS • 



«•* 




• 




•YES 


« « 


* 




IS NUHSER 


OF 


« 


♦02 • 






ITERATIONS 




9 • 


« 




TOO 8ie 




« 



*!.• 




> w «««»!t»i 



• fSTIIMTE • 



* ♦ 

• MEPME * 

• UPOATING OF • 

• MATRIX « 

• « 



IS UPOiriNG 

NUME«IC*LLV 

ST*«.E 



•¥«I.UC IN OLBF* 



,/' ■ 



•**«••*•«••*•*« 




••PIFFEIIEMCE 
*XtgMEiT$ 
• SuilEHT 



HAS FUNCTION • * 

«MU£ *t 

OECIllASEO * • 



51 
2* 



M^S^'l 



• (NO 

• IS ONE CrCLE * 
• OF ITEXJiriON ' 

• COHM.ETE • 



* 
02* 



ii^Ky* 



* CMCULilTE • • • 

• NOM OF •*•*! 00* 
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SUSRQUTIMES FHfP ANO DFNFP 



* * 
•AO • 

• m 
*•• 



• •• 

• • 

• •I 



• « 



16 



• mo 

IS GRADIENT • 
VECTOR SMALL 

eMwCH • 



**« 



«*«««** ««•!« 



DC EOUAL 1* 



FUNCTiON * 
•MMMMtCO Ml* VALUE IN ft * 
*ANO OLD OV IH* 
* OX • 



•Sr€P ARGUMENT* 
« « ALONG LINE* 
4 WITH • 

* ESTIMATED * 



•«**««**•*«•**« 



E>IT 

>*» »«♦»*♦*** *** 



•62*1 • EXIT 

• • • 
«« «**««**«4«*«««« 



GRASllNt 6 AT 



RETEAr SEAACH 



***«***•**•*•****** 



* *lfES ••• 

* 010 niECEOINC • • * 

ITERATION 'i B2* 

• FAIL TOO * i • 



***•«*•**•***** 



* SET ERROR *•**! A3 * 
*COOE EQUAL -I* * 1* 

* • *«*« 
«««*«***•«•«*•« 



* CALCULATE • 

• DIRECTIONAL * 
*OERIVATlVE or* 



«************| 



IS or ZERO 



•1» 



• I e« * 



* *Na *«•• 

• HAS FUNCTION • * > 

VALUE *l A4 * 

* DECREASED • * 3* 



* INCREASE * 

* ESTIMATED * 

* SIEPSUE • 

* * 
•***•******««** 



•J cssi IKS zh'l *C' 

M*«*«***«*M*4 



********** I* CO* 

• * 
•*« 



*U**** 



*A« *•■* RESTORE OLD ****! RO* 
» • * VALUES * * * 

*«• * • *** 

*************** 



LINE THROUGH 

X NITM 
OIRECTIOM H 
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SUaROUTINES FMFP ANO OfNFP 



A« r* It AKt 



« • 
JCOJI 



• •YES 

IS F ESIML F« •*M«**«»*M| 






«**♦•***••«*«•* I 



.a, 

« « «•• 

I* Iemch mKtmmi «s« 



• »rES 
IS f CaiML FT 



« *VES **** 

• I * * 

•I IS OT EQUM. •! U • 



• NO 



*«««*«•••«**•*• 



• • *** 

• REPLACE LIST* • * 

• S&UICH MIINT****I *4* 



• INTERMIUtTION* 



23 V 

•*4«*«««*«««««* 
«* * • 

• • • SET W • 
•*t**«l« CIMICAL • 

• • *imtei)polatii]n* 
** * polvnSmiii * 



• *YES •*•* 
IS HININIZIM6 * « 

ARGUNENI t *l 

CONPIEU * • 

• * **** 



*?.• 



•*•*«•**••••••• 



FUNCTION 

JfioKNT C «T 
FOINT 



•ea***i* oI«E 






ihmi 



• NO • • 

• I* 

36 



*****•*••* 



m 
t* CO* 
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&UaROUTlNeS fmcg aho ofhcg 



• CALCULATE • 
1* DIRECTIONAL • 

« DEKIVATXVE • 

* • 




« « 

• SAVE CURRENT* 

• AKGUMENT • 

• VECTOR ♦ 

• • 



###CAA44404t ^i^ * 



4 


VFS 


•*** 


IS 




* 


* 


01 RECTI UNM. 


»l 


An 


DERIVATIVE 


« 


* 


?* 


NEGATIVE • 




*i 


•• 




42 V 

•4« « * 

* * • SAVE NORN Uf* 
*03 •«»* CURRENT * 

• • • GRADIENT * 
•** * • 



iwTlM au 



* 

OLV* 



•••••••••••M 



• NO 
S NUMBER or * 
ITERATIONS 
TOO «IC • 



**•«••***•*«««* 



«•••••*«••«•*•• 



21 



"iJttfT' 






• VES 


• «« * 


IS 


GRADIENT • 


• * 


VECTOR SHALL •) B7* 




ENOUGH * 


* • * 




« 


*♦• i 


« 


m 


47 




« * 


*** * 




• NU 


* • * 










• • n 






*** 




•••««•••*•• **««***j 




♦ 






V 




3 


* ♦ 


I 


* 


• 






• VES 


• *• 


IS 


ITERATION • 


« * 




CYCLE 


>l A3* 


COMPLETE ♦ 


« * 




* 


**« 


• 


« 






* * 


• •• 




• NO 


* • 




«*4> •«««•** 


I* A4* 



FIK 



l« 



«•*«••*•*•«•**« 



l< OIKECTION OF* 
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SUSKOUTINES f=MC6 MtO OFNGC 



• • • tMLUC IN FT • 



*STer AKSUHENT* 



***4<(i 



*LINCUI StMCm 



It F(ll****l 89* 



COHruTE 

FUNCTION 

VALUE F «N0 

SHUIENT C AT 

nun t 



••••••««****• 



•SAVE FUNCTION* 
* VALUE IN FV * 



:«.~'aKi*«s-A». 

••• • III mm * 



„Te • 



Iwslftsi^,: 



itI*a?9«< 
c«cu 

GOHriCIE • 



****<»**««««««l 



IS or ZEU tl «3« 



•**«******l* A4 * 
20 



H » • H 


♦5 • • ^ 


IT * • 20 


• • 


* « 


• « 


•NO ••• 


•• • *VES ***• 


* «Na «**• 


HAS FUNCTION • • • 


:.5:. • iicfSS»ia[L • .t .t% 


* HAS FUNCTION « • • 


oeak!!i« ••i": 


* ViU.UE *i AA * 


* * * ENOUGH • • 1« 


* DECKEASEO • « 3* 


* ••« 


** * « *«•* 


* • «*•« 


• * 


* « 


* • 


• • ***• 


♦ * •*•• 


* • 


♦ YES ♦ ♦ 


• NO • * 


* ns 


»»»♦»»»»♦» !♦ 03 • 


—..—-,• ,j^. 


• 
• 


•««4 


••*• 


• 


42 


25 


'• 



INCKEASE 
TFFSlI 



IS STEFSIZE 

EXCEEOINCIT 

LAKSE 



•••••••••4D7***! • SET 

• • • CODE 



•ii.* 



LINE THKOUGH 

X NITH 
OIKECTION H 



14 

•*••••*«•••*••• 
**•-. i SAVE Olp • ••• 

j.4«,: vAfrii"*, :— T «: 

«„• JANO OLO OV IN* • « 

*• * 0> • •*• 

••*•••••«• 



174 



SUeftOUTINCS FHCC ANO OFMC« 



K« ft C *«C 



VM.IICS 



•NU M* 
« • • 
IS OALfA *l B^ 

NEe«TIVE • i * 





* ♦¥£$ 


:.i% * "fifssfas"' * .? «% 


* * 


• IS F £OU*L FX •••••••••••••J 


* * * FAIL TOO « • 1* 


• * 


•«« « • «•«* 


• • 



IS DX EQUAL •! A3 • 

OALFA * * 2* 



• SET EllltOR ••«•) A3 « 
•CODE EQUAL -1* • 1* 



«REPLACC OLOEft* • • 
♦*•*****•**•***>♦ S6ARCH POIW****! AS* 

♦ et CURRENT • • • 

• MINIMUM • *•• 



4 MO * • 



• « 






•• • ♦rf! 

• • « * 

•M«l IS F EOUAL F> 



***«*•«*****• I 



IS DV EQUAL 
OU.F* 






• • • KT UP * 



37 



• * ••* 
» REPLACE LAST* • • 

• SEARCH POINT****! A4* 

• BT CURRENT • i • 

• NININUM • •*• 
••*•****«**•*** 



IS NIHInlJI 

ARSUteNf 
COMPLEX 



• ARCUNENT X * 

* • 
*»* ♦«* ***♦♦♦♦♦* 



CONPUTE 

rUHCTION 

VALUE F AND 



AT 



• • 30 V 

* * *♦*•«••••••••«• 

• 'VES * * ••• 

* • CALCULATE • * • 

DOES F EXCEED *«*************|* DIRECIIUNAL O**! CO* 

FX OR FT * • DERIVATIVE • • • 

* * • OALFA • ••» 

* * «••«*•*««*«•«*• 

• • •**• 
• NO * • 
•**«******|« A3 « 

* 2» 
—»* 

38 



175 



SUiiKOUTIhie PPRCN 



* * 

♦ tNTER FPRCN • 

• * 



« • 

CAtL PERM 

* * 



VES .* IS 

,..«. THERE AN 
•. EAROR 



«*««*e3***'******« 
« * 

* COMPUTl^ • 

* ipmP2(iPi • 

* INVERSE)! * 
« * 
***************** 



♦»»*f 3**#**«#** 

« * 

* EXIT * 



•****E4********* 



•**«*«««««*«««*«>j 



( 

V 
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SUBROUTINE PERM 



***«A3** ******* 

• * 

♦ ENTER PERH • 

« * 



. OIHENSION 
♦.POSITIVE . 



.* TEST ♦ 
tPAR FOR 
PERH. 
. VECTOR 



•****c:5******«** 

« 

* ILR — 1 

« 



• ***E2**«**«***' 

[ER^l 
•*#**»»***»**«*i 



♦ CALCUL*T£ • 

* pEftHurAruN * 

* VECTOR IP2 • 

• EauivALE^r t:] * 

• n»i * 



***«E3********* 



tf ************ 



**•« *G3 ****** **^ 

• COMPOTE 

• TRANSPOSinON 
*veCTOR IPZ FOR 

• PERMUTATION 

• VECTOR IPl 

m *************** 



•***H3 ********** 

* 

* 

I€R=0, NO ERRDR* 

* 



m***-it,********* 

icR=a. N3 cRROR 
*************** 



****ji**m****** 

* * 

* EXIT * 
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SUaitOUTINES TEAS ANO OTEAS 



«•*•««*•••••••* 



«•••«*•*«*«««•« 



IS LOSS ilf • 
SIGNIFICANCE 
IN XUAUSl 



36 

•♦ • RESET • ♦♦♦ 

• * *N£U-IAUS-L AS* * * 
«*«««v«««*CO«*«|« OIHENSIUN OF*«»*» i6* 

* * •UUTPUT VECTOR* • • 
«* * X * •«* 



* * 

♦ SET NEU-N *S* 

* DIMENSION Of* 

• INPUT veCTOK* 

• X • 



•**•«*««*«•*«*« 



•«*•***•«•**«** 



«(••««««•««•««««««*«««*««« «««««**« 



♦YES ••• 
* • * 
•I C3* 

**« 



*••*••••**••*** 
« « 

• RESET Ul-U-<,* 

• T-W2t H2>U5|* 
« ANO U3«W6 * 

• « 



•* *CALCULATE Ul.* 

« • •lf2«H3iH4,HS«r* 
*A3***I* AS INITIAL A* 
^ ^ ♦RRAV-VAtUeS • 



** 



•*««««««•«**••• 



««••••«••««•««« 



ADO 1 TO I • 



•••*••«•«•«*«•* 



* SET EKRUK * 
I* PARAMETEK • 

* lER — I • 

* * 
•••*•**•***•««* 



* « 
*SET IMIN-4 AS* 

« NUMBER FOR « 

* SCRATCHED * 

* CaMf>UNENTS • 

««•********« **4 



••«*«**•«••«*«« 



******««««««««* 



c: 



«•**•***•*«•**• 

« INITIALIZE * 

• 1-5 AS * 

* COUNTER FOR * 
*Q I AGDNAL-LUOf>* 



••«*****•*««•** 

• fttSEF • 
« NEH>NEW-IHIN« 

• AS OIHENSIJN* 

• UF OUTPUT * 

• VECTOR X ♦ 
*************** 



•* ♦ SET * 

• • • lAUS*|-IHiN * 

*A6***)* AS INDEX FOR* 
« • *rHE KESULTING* 

** * vEcrtut X • 

•••*•*••*•«**•• 



IS HEUUIRED 
ACCURACY 
REACHEU 



•««•*«•««««•* ) 



*************** 

* * 

* CALCULATE • 

* M4.M^,M6 IN « 

* THE EPSILOM « 

* bl AGONAL * 
•**«*«4««««4*** 



*•****«*««*•**« 
« • 

* SPECIFY THE • 

* RESULTANT * 

* LIHIT FI'4 * 

* • 



* * 

* RESULTANT * 

* LIMIT IS ♦ 

* SPECIFIED AS* 

* FIN>XINEWI • 
**«*«««*****«** 



* •¥£$ **• 

IS LOSS UF ♦ • ♦ 
SIGNIFICANCE *t CO* 

IN W5 * • « 
» • ••« 



•**«*«*««** 



• SET ERROR • 

♦ PARAriETER ♦ 

* JEK«0 ♦ 

• * 
*«*****«•****•* 



* * 

* SET ERROR * 

* PARAMETER • 

* ItK'l « 

* * 



* CALCULATE * *•• 
« X< lAUSI * * * 

* EVENTUELL *•••} 80* 

* WITH A • • ♦ 
*SINGULAR RULE* «*• 
*•***•****«*«*« 



EXIT 
*•**«*««*••«««« 



**««•««•««**«•« 



•*«****««***«** 
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SIMDOUTINES TEUL ANO OTEUL 



CHTEO 



*• 



COMPUTE • 

•COMtl* OELT«-TERHS * 

* » • FOR • 

•• * COMPAKISON * 



• • •• 

• TES • 


••*•***•*•• 


•« 


IS MX • * 


SET EMOR 


« 


GACATER THAN «•«*««*«««*«••«}* 




• 


• * 


IER-1 


« 



•*«•* «««**•«*•« 



TKANSfOR- 

NATIQN 

FEASIBLE 



« SET EMIOII * 

• PAKAMCTEA • 

* lER— 1 • 

« * 



••*«««•«*•••••* 



«*«**«««««««««4 



•**•«••«•«••••* 



• • 

• IMITIALI2E • 

• I-L ANO n>i • 

• FOR ♦ 

• €Ut£R-iOOP • 



8 V 

• AOO 1 TQ H • 

• ANO DIVIDE * 

• LAST ♦ 
*DELTA-rERH BV* 

• 2 • 



•( •••••••••*♦•♦•»♦• 



** • • 

* * * SET J"0 AS * 
•a4*«*l* COUNTER FOR • 

* * * COMPARISUM • 
•* • WITH EPS • 



•«««««««**«*«•« 



**4t**«**4««««tt« 



*••*««*•*•«••*• 



♦85***l* AOQ I TO I * 



«««**«••*****•• 



• •¥€& ••• 
IS RELATIVt * • * 

CHANGE OF SUM «) B4* 
GREATER THAN * • • 

* EPS ♦ **• 



*N0 ♦♦♦ 

• • * 
• I CO* 

* « 



ADD I TO J « 



**«*•««*«««*«** 



EXIT 



•YES •*• 

* * * 
*l BS* 

* * « 



• m 

• SET ERROR • 

• PARAMETER * 

• IER-0 • 

• * 



EXIT 
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■Ti 



SueftOUTINES M.I *HO DALI 



•*«*•*••*•*••«* 



•*«•*••«*•*• •*• 



* PMANETER • 

• INITIM.IZE 
Ml - 



TIM.IZE 
lt?-0 



m*mm— ******** ************** ****************** ■ 



» 'Yes 

IS J EBUAL TO 



* IMiriALllE • 



•NO 
IT2 • 

THAN 



SET RESULT 
r-VALllI 



*1S. 



«AS«**)* AOO 1 TQ I 



IS J LESS 
THAN S TrESP* 



♦VES ••• 



• V6S 

* 

« 
•4 



••«•••«««••«••• 



c 



al 



jJS? 



„ • •• 

iit't Ajio *•**! ca* •* 

suiiRAcr 1 • • • •• 

man ■' - 



** *** ** * ******* 



* ************ ** ) » SET J-NOIM 



IS I EOUAt. TO *| AS* 
j-l * i * 

• * *** 



** * SET ERMM * 

* * • fMmnn » 

* * • sustIact 1 • 



« » *» « »' 



•w*'i HUiim': 



^jT^svr 



•*«•••*««**•*** 
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SimUlUTiNtS AHl ANO OAHI 



IS AHtilJI 

EQUM. la 



• • • • • COHPUTE • 
• I BO « *C0*4«|* VALtlCI 6Y • 

• 2* * • • NeANS 01= • 
**«* «« * FOWULA i * 






* COWUTE » 

• FOItNULA 2 * 



* 

QUAL TO 



• « 



• • * 
•Ol***)* 



to 



SfT ERROR 

PARAHCrBt 

IER-0 



••«•«*•*««*«*« * 



« * • 
•*** I Al • 



iU>0 2 TO I * 

•«««««4««*****« 



•««*«**«*«i*««*« 



«•*•*«*•«****•• 



• •MO 

IS NO|N |OUU ; 

• * 



•NO ••• 
« * * 
TO •• HI* 



« * 


11 


* • 




•NO 


••« 




• * 


GREATER THAN 


• 1 01* 




• * 




•«• 



• YES 



» * »» • »*♦ *»—• 











SH «f 


lULT • 




*-yMi 


1* • 




V*II2I< 


<ta * 


** 




»••«•* 



• 4 ***•*«•**««••*♦■ 



oi; 



INITlALlZe 

12-a wo 

■-0 



•res •♦•• 

•0 • 



• • Ago I ID I • 

« •*l* MtO KEPLUCC • 

• • otiii e» • 

• • OllTZ • 



• * 

• SET ERROR • 

• PARAHETER • 

• IER>I • 

• * 



**!••••« ADO 1 TO J • 



ll-«Ri<l)l «' 

J-0 » 



A&D* 



EXIT 




*« • • 

• * • tao 1 TO J • 

•M*W|* ItHO SET * 

• • « K«1«I>J ANO • 

** * •••"JiliL * 



•«i* 



•VES *«*« 



•i 
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SUMIOUTINCS «HI «H0 IMHI 



• • • ztfisiiPz • * !• ♦ • • <«vu.ni 



c 



••««••• •«•*«*«* 



•l« IfT MtUlT • 
• ir>»M.ni ; 



• • • lai-a * 

••• * * 



; 
'v. 



c 
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SUBKOUTINES ACFI «ND DACfl 



ENrCR 



•!• 



COMPUTE 
VM.|I1 (T 

• NEMS OF 

* FORMULA I • 



»*»*** *»»»>*> 00* I 



« S£T ERROR * 

• PMAMETEK * 

• IER-2 * 

• • 
*************** 



• *VES ***• 

* • * 
IS HOIM LESS •> A3 • 
THAN 1 * • 2* 



• *1iS 

IS ARCIII • 

E4UM. TO 
ARCIJI 



17 
•I «)*• 



«*•***•««•«••*« 



I* MAOIINEnlo- * 



tol * 



IS 0EI.T2 LBS 
THAN DELn 



•U-*i 



«•!••••< 



• • 
20 



•**•*••**•«• 



OELTZ-O 

* « 



* AOO 1 TO II * 



•«•*«**«**••««« 
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• • * CtUMCE SICM * 
*(0***l> OF IMASINARV* 

• • • PART • 
•• • • 



00 SIM AND * 

INVERSE 
TtaiES EXIST • 



*i 



*ai*l* EXIT 

• * * 



« (NO 

IS MAX • 

DIMENSION 
SREATJ* THAN • 
• THREE • 



• « 

*****«M*«*****I* SET IFERR-l • 



« • 

• COMrUTE _• 
•DINENSIONS OF* 
« A * 



IS raURIER 
TRAHSFOKN 
REOUIRED ^ 



* DIVIDE A SV • 

* TOTAL NUMBER* 
« OF ELEMENTS • 

* • 



* * 

* CHANCE SIGN • 
» OF INASIHARV* 

* RART OF A • 



* *VES 

IS MAX * 
01 MENS ION 
GREATER THAN • 



»>•♦«><»*** 






•••*«***•«•*••« 



ARE FOURIER 

COEFFICIENTS 

OR INVERSE 

• REOUIRED * 



• I mrmf—— ***—** 



*************** 



******** *********** 



* PUT FOURIER • 
•SUNS IN ORDER* 



** *** * ***** **** 



****«***••• 

* « 

♦ SET l«RR-l • 

* * 

• • 
*************** 



HAS INVERSE 
FOURIER 



*,l*. 
EXIT 



iJiSi***,*A ,^ 
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SUmOUTINC KHARN iWO MHMM 



•DIHimjONS I 



•SET MTHIlUn* 



• MVIDC A tV * 





( 
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SUeKOUTINES AFHf) AND OAMMt 



• SEAKCH PIVOT* 4 « 
»• OHlV using •**•! D2 • 

• PIVOT-CQLUMI* * I* 



• COMPUTE * 
I* IHTEftNAL * 

*TOLEMNCE TOL* 

• • 



•**••*«*««*•*•* 






•••«*«*•**«*«** 



•»••»••***••••* 



IS DSUN 

SHtTtO TtttN 

THe INTEIlNAt. 

TOIEMIICE TOl 



•••«•*«•< 



SIT l»UN • 
COtMl TO S • 






•A3 ••l»IU»XI!!»l„ 
* • • II AND $ 
•*• * IT IN '. 

• ••*******i 



vrii 



^U'l&i 






•wo TOflJI TO* 

• OSUM iMD • 
•STOHtTHf SUM* 

• IK TOPUr ♦ 

••••••••••••*•• 



: E»8S!riT : 

• IN Tgfiil.u** 

• TO TOMil • 

•••••••• •*•*••• 



•«•••*«••*•**•* 



• SietllflCAMCE * 



••••••*•' 



s* 



SET TOrtJI • 
(«U*t TO • 



•I *••«** «•«««««*••' 



•B7 ••I* couNfiirfaii • 
• • • Tor-Loor • 
••• * J"J-l • 

•••••••••i***** 



••*•••*•••*•••• 
*•* • • ••• 

* * * SIOU • * • 

•He ••i*coiWE$paNoiNC****l n* 

* • • COLUIIN-TE*n • • • 
••• • IN DSUN • ••• 



IS COUNTED J 
PCSITtVE 



«l 07 ♦ 
• I* 



**••*«•••••« 
UPCMTE 



• N»fMX T «NO****l M • 

• IWCX-yCCTOli^ • !• 

•••••**•••«•••• 
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SUMOUTINE! APMIl MD IMnt« 



• • • ADO 1 TO • 
*aO **l* COUNTEK ITER* 



♦NO *••• 



• StT EMUMt 



•NO 
ITER LESS * 
rMAN THE 



*«««*•««« « «•*«• 



•**••* •****•■ 



•i' 



SET CKKOIt 
•**l* p«t«NETCI> 
lER-l 



« * **** 



*««•*•«••*«•••« 



«•»«*•**•«*««*« 






««•«»«««« 



V 
EXIT 



•02 

* • 
••• 



•MO ♦••• 
• • * 

• I «2 



:ii*KiE"i &n 

• M-COIUNN • 



t-Hei 




•«•«••«•« 



««*«*Z* 



EXIT 

*************** 



• IS oniM, 

fUUAU 
KUOIEO 



imm.«yes •*•• 



*l «3 * 

***4 



*• * * 

• • • SET ISE— ISE» 
*C4***|* tm TKlNSFEtt* 

* • •¥«.«£$ IN riv* 



• SET S*»E TO • 

• PROOUtT Of 

• «EI>I '" 



■>] *N0 * 



* atlO 

• IS H Nomuic • 

• (I.E. IS tSE * 

• WSITtVEJ • 

• • 
• YES 



••«••••*•*•••*• 



•SET J>IPIV »S» • • 
I* LAB EL FOR •»«l »7» 
» PIVOT-COLUMN* * * 

* * *** 

*************** 



*************** 

* UPDATE « 

* topTipIvi as* 

* top-teIm in • 



»**»♦*♦♦<>• H i 

* SET ERROR 

* PARAMETER 

* IER-0 



*************** 
* 

• INITI 

* J*N*I 

• COUNTEL _ 

* TOP-tOOP • 
*************** 



gllZE ♦ 
AS • 
FOB • 




* MIST 
MlTH 

* BE « 



COLUMN * *TRANSFER TMIS* 

H INDEX J ***************l* COLUMN FROM * 

— * - • JUniX T TO * 

* VECTOR PIV * 

****«*••«•***•* 



CEHERATEI 



i*' 



• V6S 



••• * *YES •**« 

* * * * * « 
•07 • IS COUNTER J •! 87 • 

• • • EQUAL TO IPIV • * 2* 
*•« • * •••• 



TRANSFOR-. • 



****«**•«*«***« 



* Ci 



GEHERAT E 
NlTt At J-TH* 



* FCT 

********** 



* * 50 

* • 

* *N0 *•*• 
* MUST COLUMN * * * 

• NITH IN6c» J *l M * 

♦ »E wmmio ♦ ♦ ?• 

* * **** 



* VES 



* PREPARE * * • 

• INOEX-VECTOIl****! DO* 

♦♦»♦♦»>♦—*• »•» 



«*••••*••*••«•• 



AND SAVE 
• RESULTS IN 



*«•••*•«•♦****•• t 



IS ISE 
POSITIVE 



*•*•••**•••*••• 

* • •**• 

* COflPUTE DSUM* • * 
*USING FCT ANO«**«t BO * 
•CORRESTONOING^ • 2* 

* ROW OF T ♦ ♦♦♦♦ 
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SUeilDUTINES OTFG ANO OQTFG 



ENTFR 



*••*«****««***« 



•*••••*«**«« *• * 



««**«*•**«*«♦** 



I* INITIALIZE * 

* • 



• * * • 

♦ ♦ • AOD I TO I ♦ 
«B3***|4 AMD REPLACE * 
« * • SUfll BY SUHZ* 

*• • * 



* * 

* COMPUTf SUN2* 

* (TRAPEZniOAt* 

*ftULE, GENERAL* 

* ARGUMENT SI * 



««•«•*«**«•••«• 



• SET RESULT • 

♦ ZtI-ll«SUHl ♦ 



•••««*«******** 



♦NO *•• 
* • • 
*l 03* 



•«««• •««* «*«* « « 



IS NOIM EOOAL ••••♦♦•♦•♦♦••••|« 

TO I * ♦ 



EXIT 



*•*****«*«•**«• 
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SUSROUTINES OrfE ANO DQTFE 



* INITIALIZE • 

•»• 1-1 AHO SET • 

« HH-H/2 * 



mmmmm ••#«*•«•*« 



• • • AOO 1 TO I « 
•e3*«4>* ANO «£PLACE * 

• « « SUHl BV SUH2* 
** • • 



• • 

• COMPUTE SUH2* 

• ITRAP€2aiDAL* 



C 



*««*«**•«*««••« 



IM0IH1-SUM2* 



EIIT 



* EXIT • 
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SUflAOUTINES Q^f= AND OOSf 



*«•***•««*•«••: 



FNFFR 



*00***l* SET 1-5 



•««««#•*••••*• 



• * 

• SET HT--H/3 • 

• * 



«*«**««•**«***« 



* • • ADO 2 TO . 
*0l **t* MO KEPt4CE 

* • •SWl SY Ar — 
••• • $U»C SY 






***•**•**«**•*« 



••«*••*•«****** 



« CONPUTf SUttl* * COMPUTE AUXl* 

*****«*4**«****|* (FORHUL* 2 t*****************!* (F3RHULA 2 * 

• WITH J=3 «NI>* • WITH J-5 *flO» 

• ZllT-O} « • ZOI-SUHll * 



MtTH J«I Um* 

z<i-2i-$ymi* 



* •wo • • 

• • COMPUTE SlIMl* 

IS NOIH tFSS ♦♦••••♦•♦♦••♦♦•I* (FOflMUlA 2 • 

THAI * • ♦ MITM J-3 AND* 

* * • ZUI-OI « 



*t*f********** 



•«•«*« **«««**«4 



4********««*** 



••**••••*••••• 



« * 

♦ COHPUTf SUfl2» 

• (FOttMUlA 3 • 

♦ WITH J»* AMP* 

* ZUI-OI * 
•*•***••**«**** 



• rOHPUTE SUN2* 

* tPORHUtA 5 • 

* UITH J»4 AND* 

• 7I6)>AUXZ| * 



YES 



•NO **** 

• • « 
TO •» *0 • 

• • 2* 



• * 

• COMPUTF SU«*l* 

• (FtlRHULA 1 * 

• MITH J»2 ANO* 

• Z ( 1 I >0 1 ♦ 



**«•**«••*••*** 



SFT Z(l»=n * 



*••*«»«*«*•««*• 



****«•**** 



*** 



• CQHPUTF. 2I2J* 

• IFORtlULA 5 * 

• WITH J-2 ANO* 
« 2(A|>SUH?1 « 



• SET • 

• Z(N0IMI-AUK1* 
•ANO ZINOIH-lt* 

• >SUM7 • 

««•*•••• •• •*•• • 



• • 

• CONPUTE Utf 

• IFORHitLA ? * 

• MITH J>3 AND* 

• 2(11*01 * 

•«4****^^*«*^^« 



*••*••**«*«••** 

• * 

• COMPUTF Z(2)* 

• (FORMULA 5 * 

• HITM Ja? ANO* 

• Z(*1=SUMZI ♦ 
*•*«««*«***•«** 



•«•••••••*•••*• 

♦ SET 2(11-0, • 
•Z(3I«SUM1 ANO^ 

♦ Z(4I»$UM^ • 

♦ « 
••«««**••**•*«• 



••*•••••••**••• 



**••••••••••••• 



•«»«**«****•**« 



•****•**•**««•* 



«•*•••••• 



*<•••**•••«•*•**•*« 



^ 



12 * 

V 

•«•••*•*•*•*«•* 

FXIT 
•*««••««•***••* 



♦ * 

♦ COHPUTP 2(51* 

♦ (VOR»(UtA ? ♦ 

♦ «ITM J=5 ANO* 

♦ Z(?I«SU*11> • 
•«*«•«*••••**** 



•••• * ♦ 

•* ♦ ♦ SET * 

«« « *4^J^ 7(N0IN-1|^ 

**•* * «AUX1 • 



*( 4«««****«*****««* 



***•**««««*«•«• 



fxrr 
•••«•*«•••*•••« 



-%i< 
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SUaKOUTIHFS QSF IMO DO$F 



* « 

• • • UITH J-l+l, • 

•«• • zii-iT-sun}|« 



* * * 

IS I+l EQUAL *l 01 * 
TO NOIH * * I* 

• * 

• * *««* 

• res • • 
••*•«•••«• 1* C8 * 
• I* 

5 



C 
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SUBMMITIMES OHFG «N0 OOHFC 



INITtALIZE 
SUNZ'O 



ISii 



•NO 
JESSV 



**••««•*•«*«••« 



«• • * 

* « • AOO 1 TO I * 
•63««*l« AND KffLACE • 

• * • SUNl BY StlN2* 



•««««•««**««••« 



«**«•** *****4 ** 



lN0IH>«StlH2* 



EXIT 
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SUSROUTIflFS OHFF AW OOHFF 



CNTFR 



•««•••*«****«** 



•«*«*««*•*«*•«* 



♦NO 



•«****•••«**«** 



• SET I-l A«n ♦ 

NOIM teSS ♦•**•*»*••*♦•••»• PREPARF * 

THAN 2 • ^CONSTANTS TOH* 

* • HFRNIU RiAF* 

* « 

VFS * 



•* • • 

* • • *00 I Tfl I • 
•fi3***>» *MO REPtACF * 

* * * SUMl ev SUM7« 
•• • * 



• COMPUTF StlH2» 
•(HFRMIir hULF* 

• eOUIDISTANT ♦ 

• AftGunehT^i * 



( 





♦ Wl 


«*• 




* 


* • 


E0O*t 


TO 


*l B3* 


NDIH 


4 


* * 



• •YES 

IS NOIH EQUAL 



•*••«*•«•***•• V 



*******«**««**« 



exir 



*m**»m* ******** 
EICIT 
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SUBAQUTINeS OH-ir. *flD OQH^G 



ENTEK 

*************** 



INITIAtlZF * 

SUM2-0 ♦ 



•*****•••••***« 



LESS 
2 



**•*•«*•« 



* • 

**««r«t* IMETIALIZF • 

* I«l * 

* « 



••«*•**•«««**** 



SUMl BY SUH?* 

• ♦ 



• * 

« coMPUTF stmz* 

•IHERHITE FHR-* 
•HUL*. CENPRfL* 

• ARGUHfNTSI • 



SET RfSOUT 

zn-i)=suMi 



*««*•«*•**•««** 



Ta 



* * * 

• I B3* 

• * 



•«•**♦«*«« •«« * * 



re MniM FOUAL *«•***•••••••••)* SET RESULT * 

TO I • ^'" ♦ lfNOIM»*SU«* 

* ^ * * »••••*•*•♦•♦••• 



♦••••*«**«*4*** 
CXIT 



EXIT 
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StMROUriNFS OHSE W(0 OOHSf 



*««««*«4*i 



iUN2«0 

_ _ m 



• • 

• iCT I«| AMD * 

*COiaTA«TS FOR* 

• HnNITE RULF* 



• ADD 1 TO I * 



•********•*•*•* 



'****• 



EOyi DISTANT 

« AACUHENTSI * 



***«**««*«««*•* . 



• SET RESULT • 

* ZU-tl«SUHl * 



C 



« • 



IS NOIH EQUAL 



15!* 



•YJS 



***«*«••*** 



IN0iMI-SUH2* 



HIT 
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SU6R0UTINFS QATR AND DOATR 






** 



«***•**«****«*• 



*00 I TO J ♦ 

*ftO***>* AND Sf^T * 

♦ • «SH«SM«FCTU)*« 

•♦ • IMFH X-X*MO • 



• ♦yps 

IS 0EtT2 ♦ 
GorATCR THAN 

E 



********* 



****•**•**«•*« 

* 

COMPUTE * 

AUXm»C.5» • 

IFCTIXtl^^ ♦ 

FCTixun • 

*««•*••***•**•• 



IS J EOUAL rn 



*«« 

* * 
I 00* 

* * 

*4* 



*CI««*)* SFT lER'O *♦*•♦•♦♦**** 



**••********•** 



««*•««****•**** 



• SET H=XO-XL ♦ 



«««««««•***«*** 

• « 

* COHPUTF NEW • 

• TRAPEZOIOAL • 
*SUM AND STORE* 

♦ IT !H AUXdl* 
**•*«*********« 



****•****•*•*•*•*•* 

* 

V • 



********* 



•yps •*♦ 
* « « 

♦ » C7* 

* « 
*** 



• INITIALIZE • 

• 0«1 ANn J>1 • 

• * 



*««*«*«* 



*•««•••• *•«•««• 



I • • 10 

* ♦VES ♦*• 

* * • 

IS H EQ04L TO •» Cl« 

n • * * 

* • *•* 

♦ * 

* « 
* NU 



******••*•«**** 



• AHO 66LT2-0 • 

* * 
*«********•««** 



••*«•*•******** 
« * 

• INITIALIZE * 
«P-1, JJ>1 ANO* 

• I-l * 

• • 



***•*****«***•• 
•* « • 

* • « AOO 1 TO I * 
•AT««*I* ANO SET * 

* « * V-AUXlit. * 
•• • OCLTl»oeLT2 • 



•****••******** 



*****«**«*«**«* 



«•*«««**««***«« 



«•**«•*•«*•*•** 



•«*•*********«* 
« * 

• EXTRAPOLATE * 

• MITH F0«H,5.« 

• STORE RESULT* 
IN AUXfj) * 

-—**•*•*««* 



******••*•***«« 



«««***««******« 



*««««**«*«*«•**•< 



*«*••*«**•*••** 



•«•«***•*«*•*•* 



« *N0 ••* 

* * ♦ 
IS J EOUAL TO •! •♦• 



IS I EQUAL TO 
HOIM 



** * 

* * * 
•C7***l* SET IER-2 

* * * 



• * 
*** 



• ( •**• ««*« «•*«•*•• 




••**•***•***•** 



♦ CONPUTE TEST* 

♦ viLue oetT? • 

• • 

••«**•* *« *4«* •« 



♦ * 7 

* 

♦YES ••• 



• CONfUTf * 

• r-H«4uxrn ♦ 



INITIALIZE •*••» BO* 
SH*0 AND J-0* * • 



« NO * « 

** < *** ****>♦ CO* 

« * 
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SUBROUTINFS 0S?*OG3.. . .,QG10 AND OQG4, DOGR*. . . .IKIG}? 



fNTFrt 



• * 

• COMPOTF • 
*<XU*XLl/2 AND* 

♦ f XU- XL I * 

♦ ♦ 
»mm*»*** ******* 



« ACCUNULATE « 

• THE SUH IN • 

• FORHUtA 1 * 

• ♦ 



* • 
•HULTIPIY THIS* 

* SUH gV XI*- XL* 
•AND STORF THr* 

* RESULT IN Y • 



*•«•****•*««••• 



•*•««••«•««**« 



c 



SUBROUTINES QL 2.QL?... . .OLIO AND OOL 4,nQL8. . .. ,D0L32 






* • 

• ACCUMULATF • 

• TMF SUH IN * 

* FORMULA 1 * 

« • 



«*** * *** **** * 



•**•*«**«**•*** 



EXIT 
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suaiiouTiiiEs 0H2.aH3.....SHig am noHo.DaHif.. ...noHM 



•*•*«***•***•*« 



MCU(IW.*TC 
THC SUM IN 
KMtMUL* 1 



SUBROUTINES 0A2.013 0*10 INO Dai^.OOtR.. . . .0041? 



e«T« 



«CCUn<IL*Tf 
THC sun IN 

FORNUL* 1 



♦*«*«•**•*•*•«« 



« STORE THE » 
• RESULT IN Y ♦ 
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SUtKOUTINES 0CT3 AND ODETl 



ENTER 



»vcs 

UAL *«««*4***4C0***I* 



***«**«**•*••** 



ICR-I • 



SET EHUNt 

MIUlieTEIt 

IEI<»-J 



♦ REPLACE OYl • 



*«««*•**««««*** 



Man 1 • 



.y|S ..«. 




. w!lI«LHEJ 



* SAVE XII-ll * 

* IN A AM< • 

* Vll-ll III ■ • 






IS I fJUAl TO *»»*m»»i 



« THiOUCH l- » 



* SET RESULT * 

• IIII- * 

* ovi-ov2*ova * 

• • 



*• * * 

• ♦ • AOO 1 T0_* 
•A4<**)* STEr-COUNtet* 

« • • I * 

*m * * 



V 



c 



« SET RESULT * 

* Ovi*OVZ-OY3 • 

• * 



• . "Tg'«.!ir 



* • 
•I *«• 






««•**«••*•*•*** 



SET ERROR 

MMMETEK 

fER>4 



**•**•*•*»****««** t * 



SET RESULT • 

lah * 



******* *****••* 



-ovi*oy2*OYj» 

« * 

•4 **••****•*•* * 



*(***«**«*•**••**** 



a • 

*•«**••$••*•**• 

EXIT 

*************** 
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SU6M)UTINES 0FT3 AND XCT) 



ENTeR 



•«**«*«•*••*««* 



*•*****•* •*****)• 



SET ERROR 
PARAMETER 



«**«««****•«**« 



••*•••*•**•****•)• 



• •**• t***^^ ***« 



***«***«*•««•** I* 



PARAMETER 
IER-1 



* * 



• PREPARE • 

• OIFFEREHTIA-* 

• TION lOOP • 

« * 



* « 
» COMPUTE B BY* 

* MEANS OF * 

* FORMULA (31 * 

* • 



• SAVE • 

• V<NDtM-2» IN* 

♦ YV AND • 

* INITIALIZE * 

* i"l * 



•• 



***•*****««**«* 



*«••*•••••• 



AIH) 1 TO * 
•A6***»* STEP-COUNTER* 
* « « I AND SAVE B* 

** * IN A • 



SET ERROR 

PARAMETER 
lER-a 



• * « 
*• • 



44*******«***** 
« • 

* COMPUTE 8 BY* 

* MEANS OF • 
« FORMULA <AI * 

* « 

4 ***«•** «««••** 



*************** 

• CCMPUTE • 

• RFSULT • 

« zindimT by * 

• MEANS OF * 

• FORWLA 151 • 

« ••** ****44 *« «« 



***♦•*«*•«**•*« 



"lU'.iw : 



«*•*•*««••*«*•« 



• SET RESULT ♦ 

• 21M3IM-II • 
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* nF'^Y IN * 

« «ITXfH**1 * 

**•*****««**«** 



IS PRHT<5» 
EQUAL TO * 



•MULTIPLY H BY* 
•2, SUeTPACT !* 
« FROM IMLF« « 
* SFT TSTEP-0 * 

•*«**«««*««*«** 



•*******•*«** 



* STORE RFSIAT* 

* IM V * 



«***•«>*•«*«*«*« 



•*«***«« 



* «rPL4CF Y BY* 

* AUXtH-41 4N<1* 

* X "^Y X-'^H ♦ 

* * 

«*«*«*«**•****« 



• *YES 

IS IHLF ♦ 

GREATER THAN ♦*•**•*• 1 

10 * 



♦♦♦•♦»•*♦**♦*** 

* STORE « 

* AUXIN-?) IN ♦ 

* AUX(N-1I ANO* 

* AUXfN+^l IN « 

* AUX(N4^&1 * 
•«*•«««*****•** 



*************** 

* INTFRPOLATF * 

* V^CTlft * 

* Y(X-^H» RY ♦ 

* "(F^NS nF ♦ 

* FaRMllL* ^ * 



«*«•*«*****«*«* 



♦ PROM 

* '■n^'PUTATinN 

• OF nF«v 



IS X LAST 

POINT OF 

INTEGRATION 

» INTERVAL * 



«««««««««•*•«*« 

* STORF • 

* AUX(N-4) IN * 

* AUX1N-2) ANO* 

* A0X|N*3I IN * 

* AUX4N+5I • 
«•«**••*«****** 



*«*«*«***«****« 



**♦**♦***•*• 



••******•*«*« 
* 
CHMPOTF 
VFCTOR OF 
OFOTVATTVFS 

neoy ♦ 

««*«*4> ******* 



•(*********••*****• 



31* • 

V 

****«**««•*»«** 

FXIT 
*•*««•«****••** 



««*****•*«***•* 

* STORE • 

* AU)t(N-6l TN * 

* AUX(N-3) ANO* 

* AUX(N+U IN * 

* AUX(N+«) ♦ 
••«**«««••**««« 



*«******•*****« 

• STIR*" * 

• AIIX(N-'I IN * 

• AUX(N-3I ANO* 

• AUXIN**;* IN * 

• AtlX(N*«) ♦ 
*«*•*«*«*«***** 



* AUKIM+^I tm* 

* neiT IN X * 

* * 
«««*•***««*•*«« 



THIS IS 

LOGICAL FNO 

OF PREOICTilR- 

CORRECTOR 

PART 



V 3r»i 

*************** 

• * **«* 

• COMPUTE NFy ♦ * « 

• FRRnR-VFCrnR****! C2 • 

• AUXI 161 BY * * 2* 

• FORMULA l-l • ••*• 
••*****«*•**•** 



**♦*•**♦•♦•■ 

* 

* 

• SAVE 



IF! T 

♦ • 

*********•*♦**• 



*** 

* *♦• 

* * * 
IH ***•) 01* 

* * * 
*** 



*************** 

* • **•* 

* rn«t''UT'^ Nt^w » • • 

* rD<ir>p.yprT'>?****) n"" • 

* »UX) !ft| RV * * ?* 

* FOPMUIA 9 * •♦♦* 
*************** 
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R0NGe-«UTT4 PART Of SUBROUTINES HPCL *NT OHPCl 



♦♦♦•♦•♦•••••••••I* SAVE K IH Z • 

• * * 



* * 
Cn*<P'ITF 

VFCT1'* IF 
nC^TVATIVFS 

* OFRV ♦ 
*************** 



THIS IS 

RUNCE-KUTTA 

SUBROirriNE 
HPCL OK OHPCL 



« • 

• USF FORHIICA * 

• 11 ANO STORF* 

• RESULT AS Kl* 
« IN AUX4S} • 



?n* V 

♦•lULTiPiy n^pv* 

• BV H ANT ♦ 

♦ *;TnRF RESULT* 

♦ AS K-^ TN • 

* Atix(r) * 



*************** 

* co«PUTe • 

• FACTOR OF A • 
*IN FORMULA 12* 

* ANO STOftf rN* 

• Y • 
*************** 



*************** 

* C1«!PMTF * 

* FACna ic 4 « 

* IN F-»«M'JL4 7** 

* ANO <T1R*^ !M* 

* V * 

*«*ir *********** 



• • 

• COMPUTF * 
•ARGUMFNT X OF* 

• F ANO A IN * 

• FORMULA 12 ♦ 
*************** 



*************** 

* * 

* C'"*PUTr ♦ 

* F ANI 4 IN • 

* FPamii a ia » 



*************** 

• SFT ISM7«2 • 

• FOR RETURN • 

• FROM « 

• COMPUTATION • 

• OF OtRY * 

*****m********m 



*************** 

* * 
COHPUTF 

VECTOR ClF 
OfRIVATIVFS 

• OeRY * 
*************** 



*************** 

♦ * 

rfHP'iT*" 
■!Faiv*rivF<; 

* <1Fi»V • 
*************** 



202 V 

*************** 
♦MULTIPLY OERY* 

* ay H «Nn ♦ 

* STORE RFSUIT* 

* AS K7 IN « 

* AUX(6) * 
*************** 



7n 


*• V 






• cn^wrirr- 

• Y(X*M1 PV ■ 

• il-'SM'S "F 



• COHPUTF • 

• FACTOR Of K * 
•IN FORMULA l^* 

• ANO STGRF IN« 

• V ♦ 
*************** 



************** 

<;torf 'iFijin T* 

IN V AMP * 

2 * 

****.********** 



*************** 

• cmpurf * 

«ARGUMFNT X OF* 
* F ANO A IN • 
« FORMULA 11 * 



********* A***** 



*****-******^* 



*************** 

» SET TSII2=^ ♦ ♦♦♦ 

• FOB RFTUHN ♦ ♦ • 

• FROM •♦♦*! CO* 

• COMPUTATION » * • 

• OF DCRY * •*« 
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Mm OF suus. ma. uto ohpcl computinc oekivatives 



•«•••«•••*♦•* . .• COHPU.E MX. 

* MATRIX A * 






*••••*•***•**•• 



••*2** 



• VECTOll » BT • 

• NATKIK A ANO< 

• AOO VECTOO • 

• DEKY * 



• OERV * 



EXIT 



231 



SUMOUTIMES lavp UNO oLevr 



ENTeR 



COLUMN& 

NtTRICES B 
AND C 



«ONINI$n«TI- 

ON OF AiJOINT 

INITIAL VULUe 

mOSlEHS 



* HAS SrE(>'SIZ£ 
* CDHReCT SISN • 



• SAVE VECTOR • 
***««*«****««*4|* OF • 

•ERROR-WIGHTS* 

• IN AUXdSI • 
*************** 



*************** 



********' 



*************** 

* SET ERROR • 

* PARAMETER • 
» IHLF-12 • 

* • 

*************** 



*************** 



*************** 



****************** 



* SET ELEMENTS* 
•OF AUXTl) AND* 
*AUX<17I EQUAL* 

• TO I • 

•****«**«*«**•* 



«•«««««•«*««*«« 

• SET • 

• COLUMN-INDEX* 

• K-O AND • 
•ZERO-COUNTERS* 
•la-0 AND lC-0* 

••«•«••••«••«•• 



• INTERCHAHeE • 



• n. ANO XU • 



•*•• 

•* • 

•02 



• STEP-StZE BY* 

* lis HESiTIVE* 

♦ VALUE • 



• SMVC MUAIkf V 

• INT|ll£HiNU • 

• MATRICES B • 



•••««•*••«••«** 



»♦*•*♦ *• • **** ** 



****** ******** * 



* *ni **• 

IS X-TH ♦ • 

COLUMI) OF *i I 

MATUX C • * 
•ZERO-COLUMN* ••' 



*************** 



COMPUTE AIXUI 

Am ST-*- 

RUULt 

• NATRI.. 

******•««***••• 



fAIXU 
TORE 
LT IN 
RIX A 

•«•**•' 



• < •*•••••*«•**••••« 



**•••*••••«*•*• 

* • 

* SET ERROR • 

* PARAMETER • 

* IHIF>13 * 

* • 
**•*•**«******* 



EXIT 
**•*♦•»*• •*• ••• 



• •NO 

IS X-TH • 

COLUMN OF 
MATRIX B * 
•ZERO-COLUMN* 



***•***•***•*** 

* * 

* REPLACE * 
•AUXIl.KI BV 0* 

* AND AOO I TO* 

* COUNTER IB • 
*•««*•« •***«••• 



4i 



COMPUTE FIXUI 

ANO STORE 

RESULT IN r 

• • 

••*••••*••••••• 



*************** 

* IHITIiLIZE * 

* S«l-0.3-X« * 

* AND xtMD*XL • 



••••••**< 




ERROR-EXITS 

IN CASE OF 

HRONC INPUT 

DATA 



* *N0 

IS K'TH • 

COLUMN OF 
MATRIX C * 
•ZERO-COLUMN* 



•*****•«• 



I*****' 



* REPLACE • 

* AUXIlT.xl BY* 

* MO iOO 1 « 



•••• *********** 

• * 

• SAVE Y IN • 

• AUXI18I ANO • 

• OERV IN * 

• AUX(ll) • 
••••••»******** 



* VALUES VIXUt* 
ev FORMULA ** 



• • ••• 

• SET • • * 

• TERMINATION ***•> 01* 

* suiiCM ro * « « 

* ZERO * **• 

•*••••***«**««* 



*** 

• STORE MCA- • 
•TIVE X-pi RON* 

• OF HATUX A • 



v*n isiu 



i i. 
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SUMOUTINES LBVP AND DL6VP 



ADMINISTRATI- 
ON OF FINAL 
INITIAI. VALUE 

PROeLCH 



114 

** * SET I&U2-12 * 

* * • FUR RETURN • 
•CO***l* FKOH * 

♦ • • COMPUTATION • 
•• • OF OERy * 



««« • « 

• « * INITIALIZE • 

• Al Ml* STfP-SI2€ • 

• « • M-IXL-)Uil/l«* 

• «« • • 



110 

•* « SUBR. CELfi OR « 

* • « DGEL6 
•Sl*l GENERATES 

* * * SOLUTION OF 

«* « e*V"R IN R • 



* « 
COMPUTE 

VECTOR OF 
DERIVATIVES 

• OERY • 



• * 
«SET ISW«0 FOR* 

• RETURN FROM « 
•INITIAL VALUE* 

• PROSLEH • 

» * ♦♦♦ ♦»♦—* •** 



IS MATRIX B ••*•*••♦' 

SIN6U.AR • 



* * 
SOLVE K-TM 

ADJOINT 
INITIAL VALUE 

• PftOBLEN • 



Ul. 



> SET ERROR * 

• r«RikltETER * 

• IHIF-U • 

• * 



IS IHLF * * 

GREATER ItUN ♦♦»»•*••••••• I • EXIT 

10 * • 



«*••••«•«««««•«*•*' 



115 V 

* OUTPUT OF * 
DATA SET «, 

V« OERY. 
IHCF. NOIM 

* AND PRMT « 
*************** 



TERNINATIUN 
SNITCH ZERO 



IS IHLF EQUAL 
TO 



i»t»********»** 

• * 

• TRANSFORN * 

• VECTOR R • 

• ACCORD INS TO* 

• FORMULA a * 
** ****** * •**•*« 



• SET • 

• TERHINATIOM • 

• SalTCH TO • 

• 2ER0 • 

*************** 



*************** 

* * 
SOLVE FINAL 

INITIAL VALUE 
PROBLEH 

* * 
*************** 



• TRANSFORM • 

• MATRIX B • 

• ACCORD INC TO* 

• FORNUIA 7 • 

*************** 



• NEGATIVE • 
«ERROR-COOE OF* 

• GAUSS • 

• ELIMINATION • 

*************** 



• ^TES 

IS IHLF * 

GREATER THAN 



114 

**« 

^ • 

•I CO* 

* * 
**« 



•(**•• ♦♦•**♦***••♦• 



••*•****•**«*•* 



:-» *• '"s :*••* 


* INITIALIZE • 


:*'*• , " » §8ifs^ " . 'i 1.* 


• STEP-SIZE. • 


* SET X-XL AND* 


••• « * ***• 


• XENO-»U » 


• * 


••«•*****•****« 


* * •*• 




* VES • * 




•«•«*»•»•* 1* ■!« 




• * 




««• 




110 


* 



117 • 

V 

EXIT 

»#♦•♦*♦*•«♦•**♦ 



TOR* 



"fNlS 
VALUES 



*OEFINE 

* OF I ^ 

• VALUES 
» TIXLI-R • 

•*«««««**«««««• 



•SET I$M>1 FOR* 
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INITIAL VALUE PROaLfcM PARF OF SUBRS. L6VP ANO OLBVP 



•**•**•****«** 



***•«**««*««•*« 



400 

* SET UP N=l. • 

• *XST«X, IHLF«0* 
*♦•••••••♦•*♦♦•♦)♦ AKD • 

♦ • ER«OK-VtCIOR» 



*C0***|* AOU H TO X * 
• • • « 



SATISFACTORY 

ACCURACY 

WITHIN 10 

BISECTIONS OF 



THIS IS 

SOLUTION GF 

INITIAL VALUt 

PAOSLEHS 

GENERATED 



* STOKE V IN • 

* AUX(I) ANO * 

* OERy IN • 
« AUXU) « 



• SET ISU2-4 * 

• FOR RETURN • 

• FROM • 

• COtPUTATION • 

• OF OERY • 



•••••♦♦•«♦•*♦•• 
• * • • 

• • • • 
•01*«»l* AOO H TO X • 

* • ♦ « 
•* • * 



* SET ISHl«I • 

* FOR RETURN * 

* FKUH « 

* RUNGE-KUTTA * 

* PARI * 



• ♦ 
COHPUTE 

VECTOR UF 
DERIVATIVES 

• OERY • 



• SET ISU2»5 • 

• FOR RETURN * 

• FROM • 

• COHPUTATION • 

• OF OERY ♦ 



• COHPUTE ♦ 
VTX«H| 8Y 

MEANS OF 
RUNGE-KUTTA 

* PART ♦ 



407 V 

* SET N-2 AND • 

* SAVE Y IN • 
« AUX(21 ANO « 

* OEftV IN * 
« AUK (91 * 



* • 
COMPUTE 

VECTOR OF 
DERIVATIVES 

* OERY • 



• AOO H TU X * 

• ANO STORE * 

* V(X«ht IN « 

* AUX(2» « 
♦«**«««******«• 



• SET ISMl=3 ♦ 

• FOR RETURN • 

• FRON • 

• RUNGE-KUTfA • 

• PART « 



* SAVE Y IN * 

* AUX13) ANO * 

* OERY IN • 

* AUXIlOl • 



♦ • •AOO 1 TO IMLF* 
•85***1* ANO SUBTRACT* 

• ♦ * H FRON X ♦ 
•• • • 



•«9 «««««•««*«*« 

♦ COMPUTE * 

y(x*hT by 

MEANS UF 
RUNGE-KUTTA 

• PART • 



••*«****«*•*«*« 



***•***••••*«•• 



• ♦ 

• SAVE AUXI2I * 
*IN AUXUI ANO* 
*OIVIOE H UV 2* 

• * 



409 V 

• PERFORM TEST* 

• ON * 

• SATISFACTORY* 

• ACCURACY * 

• IFORMULA 251* 
*****«•*•«••«•* 



• SET lSHl-4 « 

• FOR RETURN • 

• FROM • 

• RUNGE-KUTTA • 

• PART * 
•**•«•«•*••««•« 



****«*«*•«*««•« 



« REDEFINE N-1* 



***•**«•*««*««* 



* * 


*•«••*««*****« « 


•YES ♦*♦ 


* COMPUTE • *•** 


IS H * « 4 


* VfX^HI 6V « « * 


SUFFICIENTLY *» 01* 


* MEANS OF *l Al < 


SMALL • « * 


* RUNGE-KUTTA ♦ • 2* 


• •♦* 


* PART • •••* 


m * 


«•••*••«««•«*** 



••••••••••♦♦*♦• 

• SET ISMl-2 • 

♦ ro« utrmn » 

* fMlH • 
« RUNGE-KUTTA • 

• PMI » 
•••••••♦*•♦♦••* 



429 



W* 



* 'NO ••• 
« * * * 

IS IHLF EflUAt •! BS* 

* TO 10 • « « 

• • •♦• 



• conruTE ♦ 
Wx*Mi a» 

MEANS OF 
RUNGE-KUTTA 

* PART • 



430 V 

* * *••• 

* SET ERROR • « ♦ 

• PARAMETER «***| Aa • 
» IMtf-lI ANO • i 3* 

• ADD H TO X • •«*• 
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INITIAL WM.UE PROBLEH PMT OF SUaRS. LBVP *N0 OLBVP 



TMIS IS START 



OF REFIWKNT 
OF START INr 

VALUCS 



RECORD OR 

HANDLE KeSULT 
VALUCS 

• • 



THIS IS START 



606 

* * • 
•00 •*■* 

* • « 



«*««***•*•••♦•« 



*A1 —i* AND AOO 



•«#**«*«••««••* 



• •NO ••*« *• « • 

IS *••••* INITIALIZE * 

t£RNIflATION •! A« • !C1*»*>* STCP-COUNTeR^ 

SWITCH ZERO • ♦ 3* • * • ISTEP-3 ♦ 

« • •*«* ♦• ♦ • 



•••********•«*• 



MEANS Of * 
A lO * 






• coHmiritliw * 

• Of OiltT * 



• IS N EdUAL TO •) 
• « • • 



600 

•res *»» 

* • * 



*4 •«••••' 



COMPUTE 
OCIIT 



• SAVE V IN 

• !u<iNI *N0 

• OEKY IN 



• INITIALIZE * 

• N-1 FOR * 

• OISPLACEHENT* 

• LOOP • 



•«•••***••*•«*• 



SAVE OCKT IN* 
AUXIllI • 



•»ES 

EQUAL TO •••««**«**C4***I« AOO I TO N • 
* ** • • 




COHfUTE • « ♦ 

--INEO » »»•••* I A*» 

CANS Of • ♦ • 

.JInULA 27 • ••• 



• HEFIN 

• MA 

• FUN 



I* AOO H TO I • 

• AND 1 TO H • 

• « 



«****«•«*«*««*•*«*< 



• CONFUTE » • . • 

• KEFINEO Y Br»»»»j A6» 

• NEANS OF • * ^« 

• FORMULA 2* • ••• 



• REPLACE • 

• AU«IM-|I ST • 
AUXNl AND • 






AUX(N««I ^ 
JkUX«N«Ti 



• AND *•* 

* * * 
IS N EQUAL TO •! C4« 

« • ••• 

* VES 

* 



• PREOlCTORlBY* 

• MEANS OF * 
« FORNULA II * 



****«« A******* • 



• STORE • 

• PREOICTOR IN* 

• AUXlUt AMD • 
•NOOIFIER n V> 



» SET ISg " 

• FOR RET. 

• FKon « 

• COMPUTATIUN • 

• OF OERY • 
••*44****««***« 



• • 
CONFUTE 

VECTOR V 
DERIVATIVES 

* OERY * 



ONPUTE • 
RECtO|_BT» 
FANS OF • 



• CONPy 

• CORREI" 

« NEAL_ _ 

• FORNULA 13 



: fL'm&Si : 



46A<>>M<> 



STARTINe 

&ALUES ARE 
NFUTED ANO 



ADO 1 TO N • 



•«*«••«••«•««•« 



CORRECTO 

FRON > 

FREOICTOR * 




IVJITI» 

OERY 



^SfS'SS 



SAVE Y IN 

iUX4H-ll A 

OfKT IN 



• * 

• STORE RESULT* 

• AS NEg • 

• eMM-VEZrOR* 

• IM UlXlUl * 




*************** 



606 

••• 
« * * * 

• AOO H TO X »•••! 00» 

• * * * 
« • *•• 



• • ** 
•COMPUTE FINAL* • 

• VECTOR Y »1 ****l 

• NEANS OF • * 

• FORMULA l« » ♦• 
*************** 
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•*4**«««****«** 



INITIAL VALUE PKOBLCH PAIIT OF SUBRS. LSVP ANO DLBVP 



THIS IS 
HALFING OF 
INCRENENT H 



• • 




628 


• • 








*•* 


SUFFi^IENTtr * 


* 
•1 


* 
Cl« 


SMAtt • 


« 


• 



* • 


601 


626 




• 




*************— 


•** A*********** 


•NO 


•«** 


•• • • 


** * « 


* 


« • 


* * • * 


* * * • 


S H TOO 


• 1 C2 • 


*C1***I*A00 t TO IHLF* 


*01***»« SET * 


SHALL * 


* 2* 


* • • • 


* • • X*0ELT-2«H • 


• 
• 


**•* 


•« « • 


•» • • 



• SCT ISII2-« • 

• F0« KETURN • 

• FRON • 

• CONPUTATIOM • 

• OF OCRV * 



601 
•YES •♦•• 



• » » 
*l A8« 



* JET ivt^-'> * 



• COMPUT. 

• OF 01 



' oil' 



CONFUTE 

VECTOR OF 

DERIVATIVES 

* OERY • 



• • 


601 


630 V 


• « 




•«*•««••««*••«« 


*M 


**«* 


* « 


ARE THERE • 


• • 


•DIVIDE H BV ?• 


ENOUCN 


• 1 C2 • 


• AM) REDEFINE* 


PRECEEOINC • 


• 2» 


• STEP-COUNTER* 


POINTS • 


***• 


•„ ISTff^,. • 



•••••«•••*••«*• 

« * 

CONFUTE 

VECTO* OF 

DERIVATIVES 

•••••••••••••V* 



• SET ISII3<1 • 

• FOR RETURN • 
•FRON HANDLING* 

• OF RESULTS • 

*************** 



IS ISXEP EVEN 



*****••*••*•••• 



•VECTOR . 

* av MEANS OF * 

• FORMULA 16 • 
*************** 



STORE T IN • 
AUR(N-2| AMD* 

8(RV IN * 
uxin»s| * 

- -«•**•***«• 



•••••••••*••*•• 

* • 
RECORD OR 

HANDLE RESULT 
VALUES 

• « 
••••***«**••*•• 



* AOO H TO H» * 

* SUBTRACT 1 • 
•FROM IHLF AND* 

* SET ISTEP-0 • 
*•*««**««****•* 



* STORE RESULT* 



*************** 



• RePLACE T BV* 

* AUXtM-41 AND* 

• X BT X-2*H * 

* * 
*«*••*•••••««*« 



IS 

TERMINATION 
SNITCH lERO 



******** 



• *Na 

IS X LAST • 
POINT OF 
INTEfiRATION • 

• INTERVAL * 



620 

*•• 

• « 

• I at* 

* • 
*** 



* I *••••*•««««**•««< 
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« bhTER OLGAH • 



.* X GT 


*. 


, hu 


♦. 10»*10 




.*.. . . 


*. 


.*' 




• . .* 






• . .* 




X 


• yes 




• * 

« * 



« * 

♦ H3 • 



**«*«U3***«****«* 



* WITH eQUATlON « 
*2» HATH. ObSCR.* 

* • 



•*«* 




•007* 




• t3 ♦.X. 




« * . 




**** 




EXITI . 




•*** 


• ♦ 


• 


♦ 6xir *x.. 


..• £3 ' 


» ♦ 


♦ 




**** 



* 


F3 •.. 


, , 










« 


• 
F3" 


X 












.« X 


GT 


'♦. 


. NO 




« 


. io*»l- 


-S) 




.♦• * • ■ 






«. 






,♦ 








*. 




,♦' 










», 


'*' 


YES 




• ♦*• 

* H3 

* 
**** 


* 
* 
* 



•••••G3****»*»*** 

. • •♦♦♦ 

♦ CUMPUTE OLNG • * ^ * 

* WITH EOUATION • X* £3 • 

•I, MATH, UfcSCR.J *»»*** 

• H3 *... 

« « . 



* S£T ERROR * 

* INDICATOR AND • 

* DLNG ♦ 

» ♦ 



241 



SUBROUTINE BESJ 



•*««***♦*«**••• 



•**«•« ««*«**4«* 



•«**4 **«••*««•* 



••««•********** 



HAS UPPEK «VES 
LIMIT FOR N **•**«••«•««• 
SEEN REACHED • 



*•*«•••«••*« **♦ 



* SET ej TO • 



** «* ••«4i*««« «** 



***•**•** ♦♦4>*«'« 



*•«*««••«•**•«• 



*••*■ 



**«**•«•*«*«•*■ 



«*«|* 85* 

« • 



••*•****«*••*•«******** 



N NEGATIVE 



r* *«•••«***•«*« 



I*** **«•*« •**4* 



EXIT 



••*«**««•«•••*••«*•••«: 



40 

*** « 

• • • 

• 62'*«*«l* 

• « • 



•••••*•••••«••« 



••**••***••«*•« 



• * 
« US£ CURRENT • 

• VALUE OF 8J • 
« FOR ANSWER « 



*«**«**4******* 



««•«••«««*«*«•* 



* INITIALIZE H* 



***•««**•**•«•* 



iSie' 



200 * 

V 
mm* *•**«««««•*«*«• 

• • * 

• C4**l* EXIT 

• • • 



••*•••*•••*•«•• 



«•***«******«** 



*•*•**«•**««*«• 



EXIT 



«*•««««•«««••*••*««•• «i 



**••*•••«*««••* 



* * 
« COHPUTE • 
« VALUES FOR • 

* ALL F • 

* * 



•CONFUTE ALPHA* 



•««*««*«**«*««4> 



•CONPUTC Nr£sr« 

• (EOUATION 41* 

* « 



• • 

• COMPUTE NEXT* 

• VALUt fOR J * 

• BESSEL * 



40 

*«« 
♦VES* • 

•»♦ Hi* 
• * * 

**« 

* 
* *«• 

• NO * * 

««««•«•«««« I* (JO* 



ESS TH 
MfeST 



REVIOUS 
VALUE 



ZOO 
**• 
•VES* * 
♦I* C4* 



242 



*******««*4|« CO* 

• * 
*** 



ibiKLUTINt t>fc*>Y 



LMtK 



**«*«*«****< •*• 



* • • :.oK Lf 1 tRMS* 



« * 

* CO* I* 

« * 



**«**««******** 



«««*•«(**«¥««** 



bXlT 
***««*«4*4**««» 



4«4*«*««***«* 



«*««««•«««««««« 



«******4***«««« 



«*««*«4«******* 

* * ««• 

*HOtTIPlV rbttH* ♦ ♦ 

• lu PRUIUCE ••*•)• HO* 

• htXl IfRH CF* * ♦ 

* tL.fb) ♦ •*• 



»*«««««*««***+*4«4**«*« 



*************** 

♦ • 

♦UHMUTfc iitt»r» 

* TtKH CF • 
*ik.o^llUfi (10)* 
« * 



*«*«:««*«««*•*** 



« • ♦ jjn Lf TIUMS* 

**« • • 



««*«**•** 4**4** 

* • *** 

* SfcT ANSWER * * • 

* EQUAL TO •*♦♦)♦ C9* 

* S£COMO TCftM « * * 

* « *«« 



♦♦«♦♦♦♦••••*••*••♦•• 



* SkT ANShER * 

* tOUAL TO *****•**•««*• 

* FIRST TERH * 

* « 



mm^m********* *4**4«*«* • 



4 • 

•K - I, YJAI =• 

• V(0», VIBI »* 

• Y(IJ * 

• * 
•4*44444*44444* 



««**«« *4**4««*« 



«9«*««««*«*i««** 



««4»«t4«««««4*» 
*4**4*««44«***« 



*»4*#4******«*«*»****» 



* • **4 

♦VfcS* * 

OCiik 16 TtKHi ♦»* CO* 



* C5****l* 

4 * * 

*44 



44444«**«*44*44 



sCOMPure • 

RECUKHENCE * 

RELATION * 

4 EO.(l) * 

*• 4444444** ««44 



* 


» 


4 


*Y 


JO 

4*v 


X 


.,ktATfck 






*)• 1 A 




THAN ^ 




* 


♦ 



««*4***4**«*4«* 

« « *44 

4KULI1PLY TbRN* ♦ * 

* TO PfUbUCC •**♦)* ft** 

* ^'fXT TLRH OF* * • 

* tg.dO) * ♦♦♦ 
4**4******44*** 



JO 

♦•♦4***4******* 

•*♦ * tllMPUIE ♦ 

• « • CUtFFIClEMS* 

• b7****]* ^OH SIN< cos* 

• * • 1EH.HS • 
•*• 4 lC-UI-(7J * 

*«*4«**««*«*44« 



*44***«44***««« 
♦CcKfUTL FIRST* 

* TtKW UF • 

* AlCUkKcNCt * 
« rLLAIlLN « 

* ^lQ.i2} ♦ 

«*»***4«« ****** 



140 
4 4 444 

• STEP INOEX 4*4 
♦FOR NEXT TERM***)* C5* 

* OF EO.dl 4*4 
4 « *«« 

*4444****444444 



444******44 4*«*«*«**««« 



* « 
*SET ANSWER TO* 
« CURRENT TERH* 
♦OF REbURRENCE* 

* RELATION * 

•«*«*444 4*44444 



4 1 4444444** 4 ***••««< 



*«4«4««4**««444 

* « «** 

* t'ifi'OTt 4*4 

* INlilAL TtKM**«*J4 (,o* 
4 Oi iCUATlcN 4 * * 

* U) • ♦*» 
**44i(:*4*4*«*«*4 



««**4*«44»«* 

tt.f>uie ♦ •♦* 

>.Cki\i} inAH * • * 

hitUKKtwCi.****!* CO* 

. t IaTUn .* * * 

tt.lij * •♦* 

:»*«4«*«4**** 



*«* *««»4>4«4 4444444 
4*4 

* C9**i* EXIT 

* 4 4 

44* ******* 4* ***4«* 



243 



^UdKUUTlNt LXPl 



•*«««*«««««*« 



**««**«*i*«4**** 



I^ J( CUUAl 1C 



••««*4*«** 



««««««*«*«*•««« 



* * 

* LVALUAU * 
« »%LVNOHiAL * 

• ACfRaXI#tATIljN>^ 

• » 



•V£S 
IS X CRcATEft ♦♦•••*♦••••*• 
THAN 4 * 



• * 
*KL:^ULT hCUALS« 
« HCjLYNLMIAL * 

• VALUE PLUS * 
^ALUGfAtfSCXll* 







• 


« 


*N0 


30 
• • 


*i 


X 


LESS 




•1* 60* 


THAN 


-4 




* 


* * 



cXlT 



«««««««««*«««*««««**««« 



**«***«4* **«**• 



••«*««««»«* 



**««•«**•«*«*** 



4***«***4**««** 



£XIT 



fcXIl 



«4«**«4 ******** 



* SET ARGUMENT* 
«EQUAL TO 4./X* 

• • 
«**«*«****4«*** 



« • 

* eVALUATt: * 

* POLYNOH]AL • 
*APPftaxiHATIlJM* 

« * 



•«**««4**4*4 •« « 

•RESULT EfiUALS* 

* POLYNOMIAL • 
« VALUE TIMES • 

♦ WP<-XI • 



«••******•«*•** 



***• •«*••••«•*« 



246 



ibttkuUl iNr. blLl 



*»♦♦****■♦*****■* 



«*««***««4««« 



«*«««««**«**««* 



•««*•*«««**««** 



* * 

» id ii t^UAL* 

* Y .;lu<il to ♦ 

*■ CLjIaI ♦ 



♦ ♦ 

* tVALUATt •* 

* PlLYKC<M1AL * 

♦ FCK SI * 



« * 

t iVALljitl;- * 
* PCLYhuM aL • 



««*•**«******«« 



****«*«■ 



♦ * 

♦ tVALJATt • 

♦ APfKOXlMATltJN* 

♦ Foi» CI ♦ 

♦ ♦ 



EXIT 

*«*****« ******* 



* * 

* * 



♦*♦*♦*♦♦*♦♦♦•*• 

cxir 

««44«^*«*«***«* 



LALCUlAIL: 
t.Y*Vl 



***** 



4««4 J9«** 



« VflLOLS KJK * 
* hcLAriVL * 

« AlitiUM^nl X * 



lXII 



247 



:>LflKLOI INL CS 



tMTtK 



***«*««**««*4** 



«**«****«*«««** 



* « 

♦VtS • itl t fctUAL ♦ 

* ♦ s b*iuAL ra ♦ 

» SiN(Z) ♦ 



«****«***•**»*« 



************* 



* « 

* fcVALOAId • 

* POLYNUhlAL ♦ 

* l-Lh A AM, 6 * 

*************** 



*************** 

♦ * 

♦ St I 5 LCIUaL * 

♦ TL l*C * 

♦ RfcPLAtt i BY* 

♦ 2*2 ♦ 
**********4**** 



*4.**i********V* 



*v.4********t*A* 



**************** 

* * 
« EVALUATE * 
*APPftOX]HAT]UN* 

*FOR c am; for* 

* S ♦ 
*************** 



*************** 

* * 
« LALCULATL ♦ 

* C=,5C ♦ 

* /*(C*AtS*bl * 

* « 
******* ******t* 



*********** 

LXJT 
*************** 



t********4***t* 

* * 

* CALCULATt * 

* i=.jt * 

* i*(ii»A-i*cl » 

* » 
*>(.************* 



ti^************* 

Lxn 



248 



SU6RUUT1HES CtL 1 ANU DCfcLl 



ENTER 



***«*«***•*«•*• 



«*4«4*«**««««*«)*KESULT EQUALS* 

* INdNlTV * 

♦ * 



•«*««*«4******* 

• ♦ 

• INITIALIZE « 

• ARITHHETIC • 
•AMD GEOMETRIC* 
*NCANS ARIfStS! 



EXIT 



« * 

• SET AARI-ARI* 

"OMPUIE * 

-ARI+6EQ • 

* * 



•♦**♦•*♦*•*•*•• 



**•****««**•*•* 



* ♦ 

• COMPUTE NEW • 

* GEOMETRIC • 
« MEAN GEO • 

♦ • 
44 ••*•** ****««« 



••»«4«44*4***** 

EXIT 
«*••••* •«»4**** 



**«**44*«4**«** 

• 4 44* 

• 4 • • 

4 REPLACE ARI ••*4| A34 

• ev ARI/2 • %** 

••44444**»«*#** 



249 



:iUttROUTlNES CfcL2 *NJ OCEL^ 



ENTER 



* 

IS CK EQUAL ***«•**•««•*«] 



*«***««*****«*« 



*«4>«*«*«*«««««« 



* iNrriALUE * 

• AKITHHETIC • 
•ANO CEOHeritlC* 
•HEAMS ARI,GEO« 



♦YES * * 

• ♦ ♦ 

EOUAL TO ••**♦*♦♦*♦•*••* i^KtSUlT tUUALS* 

O * * A « 

• • ♦ 



*( «*«*******•«*«**« 



«•** * «** • ***•«* 



**•••••***«• *•« 



«««*•«*«««*«««• 



*«*««*«**««««** 



EXIT 



«• • • 

* * • COMPUTE • 
«A4«**)«W-2IU+AA*GeOI* 

* • • SET AA-AN * 
•« « • 



EXIT 



***4* *«*•««••«* 

• SET AAKl'AKI* 

• COMPUTF • 

• AftI>AftI«GEO « 

« • 



•***•***«***•«« 



*♦*****•******« 



*•••*««•••<•«*•« 



•*••*•*«******* 



**««*«4-«*«**««* 



• COHPUTE NEW ****! A4* 

« Geo » « « 

* * ••« 



********•**«« 



♦•«*•*«**•**# 
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sueouuriNts ttii (nu DtLU 



ENTEK 



•YES 



««*««*•«*•«**** 



• IS 4RSUMEN1 X »»»«»•»«••«••«•»» RESULT li J »•«»*••••••♦»•»« I » 6X11 

• EUU»t ^ • J J * »»«««...*..... 

• NO 



* •YES ♦ * 

, IS • •S€I UP BESUH* 

CIMPLCHFNTAIIY **»«««*M«>****I* VUUE «S • 
• iSobuLUS » • ELE«ENI*R» • 



• JiIgEO^ANO I »**• • EVAlUAIt • 

• VALUES • • • • COKttECr SIGN* 

• ANGLE. PIN • •• l.,.»„»»»...,I 



»"♦ • COHPUTE KEN • • ^,,, 

•A4***l* Adl PREPARE * • EXIT 

%.♦ : "* "" «" : 



•**«♦***•**«*«* 



VALUE} 



* * ■ 1 

* NO * 



« • ♦ * *** 

• UPDATE * * * • ♦ 

• MALUe PIM • • VAtOE • • • 

• • * * •** 
««tttt**«*««*«**4 *************** 



« • *«« 

• * • • 

• COftPUTC HEV ••••I A4* 

• ceo ♦ « * 

• 4 *«« 
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SUBRUUIINfeS ELI2 Af«D UbLI2 



FNTFR 



*«**•*«***•*««* 



«*««* *««•«**««« 



INITIALIZE « * 

VALUES C.O * * exiT 



* VALUES C.O * 
« FOR SUH OF * 



• *irES • ♦ 

IS * *S£T UP kESULT* 

COHPLbHENTARY ««««****«««*««« J*AS ELtHbNTAKV* 

HUOULUS O * * INTbGAAL * 

* * • ♦ 



7 V 4 V 

* INITIALIZE « ** • * 

* AA.AN.K AND 

* ARdUHENT 
VALUtzS 



« ANG,PIN,ISI * 



«*«•**«•**««*«* 



* INiriALIZE « * 

♦ AHirHHETIC ♦ • tXIT 
•AND GEOMETRIC* * 

*HEANS ARI.^a* ***********0*** 



«««**««« •«•«««* 



• COMPUTE NEb * 

* RvSGEO AND • 

• NEW ARI * 

* * 



• « 

* UPDATE • 

* ARGUMENT * 
«VALUES SUH OF* 

• SINES ANJ * 



*****««««*««*** 



* VALUfc ♦ • ♦ 

• « *** 



* UPDATF * •** 

* ARGUHENT * « * 

* VALUES AND ****> UO* 

* CaMPUTF NEW * • * 

* GEO * ««* 
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SUSROUTINES JECF AND OJELF 



ENTER 



»YES • • 
• * EVALUATE * 
IS NQOULUS **«*****«******|* SN>CN«DN AS •«*«««*««•*•••«« I* 
EOUAl. I • ' --..I— 



***•«•«•••«•««« 



• ELEHENTAIty • 

• FUMCTIONS • 



«««**«««««*•««• 



*VES * PERFIMN « 

__ • „ • JACOBIS » 

IS MOOUiUS ••****•*»•*••« I* MODULUS * 

CHEATER 1 • • TRANSFOR- • 

* • HATION * 

• * • 

* NO • 



• (•••••••••*•«•****•••*«**••••**•* 



««*«««*< 



• • •CALCULATE AM)* 
•A««**I*SIDRE NEW ARI* 

• * * AND CEO • 



IS ARI-GEO * • 

SUFflCIENTlY •••*•«***•••••• I* 

SMALL • • 



* «** 

♦ NO ♦ * 

*«**««**««|« A4* 



•*••*•««••••«•« 



••••*««*) 



INTERCHANGE • 
CM AND OH • 



••«***•••«***«• 



* SCHEHe FOR • 

• ANCLE « 



*************** 

* * 

* CALCULATE • 

* SN.CM FROM * 

* COfANCENT « 

* ANI>LITUOINIS* 
*************** 



«| ****************, 

* 
V 



*************** 
ERIT 



*N0 **************^ 

* * 

S *************}* EXIT 

• « 

* *************** 

*** 
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